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Abstract 

Extensive experimental data from high-energy nucleus-nucleus collisions were recorded 
using the PHENIX detector at the Relativistic Heavy Ion Collider (RHIC). The 
comprehensive set of measurements from the first three years of RHIC operation 
includes charged particle multiplicities, transverse energy, yield ratios and spectra 
of identified hadrons in a wide range of transverse momenta (pt), elliptic flow, 
two-particle correlations, non-statistical fluctuations, and suppression of particle 
production at high px- The results are examined with an emphasis on implications 
for the formation of a new state of dense matter. We find that the state of matter 
created at RHIC cannot be described in terms of ordinary color neutral hadrons. 
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1 INTRODUCTION 



1.1 Historical Introduction 

A recurring theme in the history of physics is the desire to study matter under 
extreme conditions. The latter half of the twentieth century saw this quest ex- 
tended from 'ordinary' atomic systems to those composed of nuclear matter. 
Even prior to the identification of Quantum Chromodynamics (QCD) as the 
underlying theory of the strong interaction, there was considerable interest in 
the fate of nuclear matter when subjected to density and temperature ex- 
tremes [1,2,3]. Particularly intriguing was the suggestion that new phases of 
nuclear matter could be associated with a corresponding change in the struc- 
ture of the vacuum [4]. These considerations gained additional impetus with 
the realizations that a) QCD was the correct theory of the strong interaction, 
b) the phenomena of quark confinement was a consequence of the nonpertur- 
bative structure of the vacuum and c) this vacuum structure is modified at 
high temperatures and/or densities, suggesting that quarks and gluons under 
such conditions would be deconfined. Taken together, these facts suggest that 
QCD is a fundamental theory of nature containing a phase transition that is 
accessible to experimental investigation. 

It is quite remarkable that this understanding was achieved very early in the 
development of QCD. Collins and Perry noted in 1975 [5] that the reduction 
of the coupling constant at small distances indicated that the dense nuclear 
matter at the center of neutron stars would consist of deconfined quarks and 
gluons 2 . Their treatment focused on the high-density, low-temperature regime 
of QCD, but they did note that similar arguments might apply to the high 
temperatures present in the early universe. An extensive review by Shuryak 
in 1980 [7] is the first to have examined the high-temperature phase in detail, 
and is also notable for proposing the phrase "quark-gluon plasma" (QGP) to 
describe the deconfined state: 

When the energy density e exceeds some typical hadronic value (~ 1 GeV/fm 3 ), 
matter no longer consists of separate hadrons (protons, neutrons, etc.), but 

Email address: PHENIXSpokesperson: zajc@nevis.columbia.edu 
(W.A. Zajc). 

1 Deceased 

2 In fact, prior to the development of QCD the quark hypothesis raised serious 
issues concerning the stability of neutron stars [6]. 
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as their fundamental constituents, quarks and gluons. Because of the ap- 
parent analogy with similar phenomena in atomic physics we may call this 
phase of matter the QCD (or quark-gluon) plasma. 

Developing a quantitative understanding of the deconfining phase transition 
in hadronic matter and of QGP properties has proven to be a challenging task. 
While simple dimensional arguments suffice to identify both the critical energy 
density Ec ~ 1 GeV/fm 3 and the associated critical temperature T c ~ 170 
MeV, these values also imply that the transition occurs in a regime where the 
coupling constant is of order unity, thereby making perturbative descriptions 
highly suspect. 



Progress in understanding QCD in the extremely non-perturbative domain 
near the critical temperature has relied on an essential contribution by Creutz 
[8], who showed that numerical implementations of Wilson's lattice formu- 
lation [9] could be used to study phase transition phenomena. This work, 
together with the continued exponential increases in computing power, stim- 
ulated the development of lattice QCD, which in turn has led to detailed 
investigations of the thermodynamic properties of quarks and gluons [10]. 




i ■ r 



3 flavour 
2+1 flavour 
2 flavour 



2.0 



2.5 



3.0 



T/T r 



3.5 4.0 



Fig. 1. Lattice QCD results [11] for the energy density / T 4 as a function of the 
temperature scaled by the critical temperature Tc- Note the arrows on the right 
side indicating the values for the Stefan-Boltzmann limit. 

Lattice QCD predicts a phase transformation to a quark-gluon plasma at a 
temperature of approximately T « 170 MeV ~ 10 12 K, as shown in Fig. 1 [11]. 
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This transition temperature corresponds to an energy density e ~ lGeV/fm 3 , 
nearly an order of magnitude larger than that of normal nuclear matter. As 
noted above, this value is plausible based on dimensional grounds, since such 
densities correspond to the total overlap of several (light) hadrons within a 
typical hadron volume of 1-3 fm 3 . No plausible mechanism exists under which 
hadrons could retain their in vacuo properties under these conditions. Lattice 
calculations also indicate that this significant change in the behavior of the 
system occurs over a small range in temperature (~20 MeV), and suggest that 
the change of phase includes the restoration of approximate chiral symmetry 
resulting from greatly reduced or vanishing quark constituent masses. 



T 



150 MeV 



QGP 




m Uid = 0; m s = oo 






_ Tricritical Point 




Hadronic 








<\ , 


\^ 2SC 



1 GeV 



Nuclear Matter ■ 



fj, 



Fig. 2. Theoretical phase diagram of nuclear matter for two massless quarks as a 
function of temperature T and baryon chemical potential [i [12]. 



In the limit of massless noninteracting particles, each bosonic degree of free- 
dom contributes §qT 4 to the energy density; each fermionic degree of freedom 
contributes | this value. The corresponding "Stefan-Boltzmann" limits of the 
energy density Ssb for the case of 2(3) active flavor quark-gluon plasma is 
then 



£sb 



7 7T 2 7T 2 

{2 f ■ 2 S ■ 2 q ■ 3 C - + 2 S • 8J— T 4 = 37 — T 4 

8 J 30 30 

7 7T 2 7T 2 

{3 f • 2 S ■ 2, ■ 3 C - + 2 S ■ 8J— T 4 = 47.5 — T 4 

8 J 30 30 



(1) 
(2) 



after summing over the appropriate flavor, spin, quark/ ant iquark and color 
factors for quarks and spin times color factors for gluons. The large numerical 
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coefficients (37 and 47.5) stand in stark contrast to the value of ~3 expected 
for a hadron gas with temperature T < Tq, in which case the degrees of 
freedom are dominated by the three pion species 7r~, 7r°, 7r + . 

The exact order of this phase transition is not known. In a pure gauge the- 
ory containing only gluons the transition appears to be first order. However, 
inclusion of two light quarks (up and down) or three light quarks (adding the 
strange quark) can change the transition from first order to second order to 
a smooth crossover. These results are obtained at zero net baryon density; 
dramatic changes in the nature of the transition and in the medium itself are 
expected when the net baryon density becomes significant. A schematic version 
of the phase diagram for an idealized form of nuclear matter with vanishing 
light quark (up and down) masses and infinite strange quark mass is presented 
in Fig. 2 [12]. For sufficiently large values of the baryon chemical potential \x 
this system exhibits a first order phase transition between hadronic matter 
and QGP, along with a tricritical point below which the transition becomes 
second order. However, non-zero values of the light quark masses dramatically 
alter this simple picture: The second order phase transition denoted by the 
dashed line in Fig. 2 becomes a smooth crossover, and the tricritical point 
correspondingly becomes a critical point designating the end of the first order 
transition found at higher values of fi. For example, recent calculations [13,14] 
indicate that the transition is a crossover for values of /i <~ 400 MeV. Given 
that both theoretical arguments and experimental data suggest that nucleus- 
nucleus collisions at RHIC (at least near mid-rapidity) are characterized by 
low net baryon density, we will restrict our attention to this regime, while 
noting that the predicted smooth nature of the transition in this region in- 
creases the experimental challenges of unambiguously establishing that such a 
transition has occurred. We also note that while Fig. 2 shows that the region 
of low temperature and high baryon density is expected to show a transition 
to a color superconducting phase of matter, this regime is not accessible to 
RHIC collisions and will not be discussed further. 

While the lattice results plotted in Fig. 1 show that the energy density reaches 
a significant fraction (~ 0.8) of the Stefan-Boltzmann values in the deconfined 
phase, the deviation from Esb, and the reason for the persistence of that 
deviation to the highest studied values of T/T c , are of great interest. For 
instance, Greiner has noted [15] that "in order to allow for simple calculations 
the QGP is usually described as a free gas consisting of quarks and gluons. 
This is theoretically not well founded at T m T". In fact, analysis of the 
gluon propagator in a thermal system [16,17] has demonstrated that effective 
masses of order g(T)T are generated, suggesting that the relevant degrees of 
freedom are in fact massive near Tq- m g ph T c could be generated by gluons. 
Especially interesting is recent work which indicates that both heavy [18,19,20] 
and light [21] flavor states may remain bound above Tc, calling into question 
the naive interpretation of e(T) as an indicator of the explicit appearance of 
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quark and gluon degrees of freedom. This is supported by explicit calculations 
of the spectrum of bound states above Tq [22] which predict a rich structure 
of states that belies a description as a weakly interacting parton gas. 

To emphasize this point, consider the standard measure of the degree of cou- 
pling in a classical plasma, obtained by comparing the relative magnitudes of 
the average kinetic and potential energies: 

r ^ <^>. (3) 

(Ekin) 



In the case of the QCD plasma the mean inter-particle spacing should scale 
as some numerical coefficient times 1/T. Naively, this gives a mean potential 
energy (V(r)) ~ a s (T)(l/r) ~ a s (T)T, leading to 

a s (T)T 

T ^ a,(T). (4) 



Any reasonable estimate for the numerical coefficients leads to T > 1, which is 
the condition for a "strongly-coupled" plasma. In reality, the screening present 
at such densities (or equivalently, the generation of effective gluon masses) 
modifies the mean potential energy to (V(r)) ~ g(T)T, which only increases 
the estimated value of T [23]. Considerations such as these have led some 
authors [24,25] to denote quark-gluon plasma in this regime as "sQGP" for 
"strongly interacting QGP" . 

It is worth noting that this state of affairs has been anticipated by many 
authors. Whether the argument was based on the divergence of perturbative 
expansions [26], on phenomenological descriptions of confinement [27], on the 
development of effective gluon masses from plasmon modes [28] or on general 
principles [15], it is clear that the QGP near T c should not be regarded as an 
ideal gas of quarks and gluons. 

How high a temperature is needed not just to form a quark-gluon plasma, but 
to approach this "weakly" interacting plasma? A calculation of the pressure 
of hot matter within perturbative QCD [29] is shown in Fig. 3. The pressure 
result oscillates significantly as one considers contributions of different orders. 
These oscillations are an indication that the expansion is not yielding reliable 
results. However at temperatures approaching 1000 times of T c (~ A^g), they 
appear to be converging toward the Stefan-Boltzmann limit (asymptotically 
free partons). It is interesting that in considering the highest-order term, the 
results are still nonconvergent though one seems to approach the lattice cal- 
culated pressure. Unlike the case of single parton-parton scattering at zero 
temperature, the infrared problems of finite-temperature field theory prevent 
further analytic progress even for very small values of the coupling constant 
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Fig. 3. Perturbative QCD results for the pressure as a function of temperature at 
various orders normalized to the Stefan-Boltzmann value ps b [29] . 

[29,30,31]. 

The goal of relativistic heavy ion physics is the experimental study of the 
nature of QCD matter under conditions of extreme temperature. A great em- 
phasis has been placed on "the discovery of the quark-gluon plasma", where 
the terminology "quark-gluon plasma" is used as a generic descriptor for a 
system in which the degrees of freedom are no longer the color neutral hadron 
states observed as isolated particles and resonances. This definition is lim- 
ited since high-energy proton-proton reactions cannot be described purely in 
terms of color-neutral hadrons, but rather require analysis of the underlying 
partonic interactions. The hoped-for essential difference in heavy ion collisions 
is the dominance of the partonic-level description for essentially all momentum 
scales and over nuclear size distances. Beyond this simple criterion, in order to 
characterize the produced system as a state of matter it is necessary to estab- 
lish that these non-hadronic degrees of freedom form a statistical ensemble, 
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so that concepts such as temperature, chemical potential and flow velocity 
apply and the system can be characterized by an experimentally determined 
equation of state. Additionally, experiments eventually should be able to de- 
termine the physical characteristics of the transition, for example the critical 
temperature, the order of the phase transition, and the speed of sound along 
with the nature of the underlying quasi-particles. While at (currently unob- 
tainable) very high temperatures T ^> T c the quark-gluon plasma may act as 
a weakly interacting gas of quarks and gluons, in the transition region near T c 
the fundamental degrees of freedom may be considerably more complex. It is 
therefore appropriate to argue that the quark-gluon plasma must be defined in 
terms of its unique properties at a given temperature. To date the definition is 
provided by lattice QCD calculations. Ultimately we would expect to validate 
this by characterizing the quark-gluon plasma in terms of its experimentally 
observed properties. However, the real discoveries will be of the fascinating 
properties of high temperature nuclear matter, and not the naming of that 
matter. 



1.2 Experimental Program 

The theoretical discussion of the nature of hadronic matter at extreme densi- 
ties has been greatly stimulated by the realization that such conditions could 
be studied via relativistic heavy ion collisions [32]. Early investigations at the 
Berkeley Bevalac (c. 1975-1985), the BNL AGS (c. 1987-1995) and the CERN 
SPS (c. 1987-present) have reached their culmination with the commissioning 
of BNL's Relativistic Heavy Ion Collider (RHIC), a dedicated facility for the 
study of nuclear collisions at ultra- relativistic energies [33]. 

The primary goal of RHIC is the experimental study of the QCD phase tran- 
sition. The 2002 Long-Range Plan for Nuclear Science [34] clearly enunciates 
this objective: 

. . . the completion of RHIC at Brookhaven has ushered in a new era. Studies 
are now possible of the most basic interactions predicted by QCD in bulk nu- 
clear matter at temperatures and densities great enough to excite the expected 
phase transition to a quark-gluon plasma. As the RHIC program matures, 
experiments will provide a unique window into the hot QCD vacuum, with 
opportunities for fundamental advances in the understanding of quark con- 
finement, chiral symmetry breaking, and, very possibly, new and unexpected 
phenomena in the realm of nuclear matter at the highest densities. 

The RHIC accelerator and its four experiments were commissioned and brought 
online in the summer of 2000. The initial operation of both RHIC and the 
experiments has been remarkably successful. In these first three years the ac- 
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celerator has collided, and the experiments have acquired data on, Au+Au 
collisions at five energies, an essential p + p baseline data set, and a critical 
d+Au comparison. The analyses of these various systems have resulted in a 
correspondingly rich abundance of results, with over 90 publications in the 
refereed literature. 

It is therefore appropriate to reflect on the physics accomplishments to date, 
with a particular emphasis on their implications for the discovery of a new 
state of matter. At the same time, it is essential to identify those features of 
the data (if any) that are at odds with canonical descriptions of the produced 
matter, to specify those crucial measurements which remain to be made, and 
to outline a program for continued exploration and characterization of strongly 
interacting matter at RHIC. The PHENIX collaboration [35] has performed 
such an assessment; this document represents a summary of its findings. 

The PHENIX Conceptual Design Report [36], submitted to BNL/RHIC man- 
agement on January 29th, 1993, outlined a comprehensive physics program 
focused on the search for and characterization of new states of nuclear matter. 
The measurement of electromagnetic probes and high-transverse-momentum 
phenomena formed a major thrust of the proposed program. It was also re- 
alized that the measurement of global variables and soft identified hadron 
spectra in the same apparatus was essential to the goal of understanding the 
evolution of the produced matter over all relevant timescales. These diverse 
criteria required combining an unprecedented number of subsystems together 
with a high-bandwidth trigger and data-acquisition system into an integrated 
detector design. Particular attention was given to minimizing the conflicting 
design criteria of the central arm spectrometers, with their requirement for 
minimal mass in the aperture, and those of the muon spectrometers which 
require maximal absorption of the incident hadron flux. The data acquisition 
and trigger system was designed to accommodate the great variety of inter- 
action rates and event sizes provided by RHIC. Every effort was made to 
provide for future upgrades, both in the geometry of the experiment and in 
the architecture and design parameters of the read-out system. 

The published PHENIX results of Au+Au collision at a center-of-mass energy 
per nucleon pair, v /s^7, of 130 GeV [37,38,39,40,41,42,43,44,45,46,47,48] and 
at v /5^7 = 200 GeV [49,50,51,52,53,54,55,56,57,58,59], p + p collisions at 
y/s = 200 GeV [60,61,62,63], and d+Au at = 200 GeV [64,65] clearly 

demonstrate that PHENIX's goal to make high-quality measurements in both 
hadronic and leptonic channels for collisions ranging from p + p to Au+Au 
has been realized. A summary of these results illustrates this point: 

• First measurement of the dependence of the charged particle pseudo-rapidity 
density [37] and the transverse energy [38] on the number of participants in 
Au+Au collisions at ^s^n = 130 GeV; systematic study of the centrality 
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and ^snn dependence of dEx/drj and dN c h/drj [59]. 

Discovery of suppressed production for 7r°'s and charged particles at high 
p T in Au+Au collisions at y / s/v r /v = 130 GeV [39] and a systematic study of 
the scaling properties of the suppression [47]; extension of these results to 
much higher transverse momenta in Au+Au collisions at y / s/v r /v = 200 GeV 
[49,53]. 

Co-discovery (together with BRAHMS [66], PHOBOS [67] and STAR [68]) 
of absence of high-p T suppression in <i+Au collisions at y/s NN = 200 GeV 
[64]. 

Discovery of the anomalously large proton and anti-proton yields at inter- 
mediate transverse momentum in Au+Au collisions at ^/snn = 130 GeV 
through the systematic study of 7r =t , K ± , p and p spectra [40]; study of the 
scaling properties of the proton and anti-proton yields in Au+Au collisions 
at ^5^7 = 200 GeV [52]; 

Measurement of A's and A's in Au+Au collisions at y's/wv = 130 GeV [43]; 
measurement of 0's at ^s^n = 200 GeV [69] ; measurement of deuteron and 
anti-deuteron spectra at y/s NN = 200 GeV [58]. 

Measurement of Hanbury-Brown-Twiss (HBT) correlations in 7r + 7r + and 
it~it~ pairs in Au+Au collisions at y/s^N = 130 GeV [41] and 200 GeV 
[56], establishing that the "HBT puzzle" of i? ou t ~ R s i^ e extends to high 
pair momentum. 

First measurement of single electron spectra in Au+Au collisions at ^/snn — 
130 GeV, suggesting that charm production scales with the number of binary 
collisions [42]; measurement of centralty dependence of charm production 
in Au+Au collisions at y/sj^ = 200 GeV [57]. 

Sensitive measures of charge fluctuations [44] and fluctuations in mean px 
and transverse energy per particle [45,55] in Au+Au collisions at v /saT/v = 
130 GeV and 200 GeV. 

Measurements of elliptic flow for charged particles from Au+Au collisions 

at v /5^7 = 130 GeV [46] and 62 GeV to 200 GeV [70] and identified charged 

hadrons from Au+Au collisions at y/s^N = 200 GeV [50] . 

Extensive study of hydrodynamic flow, particle yields, ratios and spectra 

from Au+Au collisions at y/s]^ = 130 GeV [48] and 200 GeV [54]. 

First observation of J ftp production in Au+Au collisions at y/sj^N = 200 

GeV [51]. 

Measurement of the nuclear modification factor for hadrons at forward and 
backward rapidities in g?+Au collisions at y/sNN = 200 GeV [65]. 
First measurement of the jet structure of baryon excess in Au+Au collisions 
at v /5^7 = 200 GeV [71]. 

First measurement of elliptic flow of single electrons from charm decay in 
Au+Au collisions at y/s]^ = 200 GeV [72]. 

First measurement of direct photons in Au+Au collisions at v /sJv r /v = 200 
GeV [73]. 

Measurement of crucial baseline data on 7r° spectra [60], direct photon pro- 
ducion [63], and J/ip production [61] in p + p collisions at y/s = 200 GeV. 
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• First measurement of the double longitudinal spin asymmetry All in tt° 
production for polarized p + p collisions at y/s = 200 GeV [62]. 

These publications encompass physics from the barn to the picobarn level; 
their very breadth precludes a detailed presentation here. These data, together 
with a rich program of future RHIC measurements, will allow us to address 
many of the features that would characterize a quark-gluon plasma: 

• Temperature 

• Parton number density 

• Energy density 

• Opacity 

• Collective behavior 

• Thermalization leading to the quark-gluon phase 

• Deconfinement 

• Number and nature of degrees of freedom 

• Recombination of quarks and gluons to form final-state hadrons 

• Chiral symmetry restoration 

• Time evolution of system parameters 

• Equation of state 

• Color and thermal transport properties 

• Critical behavior 

As emphasized above, the present PHENIX data set from RHIC runs in year 
2000 to 2003 already provides an extensive set of measurements on global 
variables: (transverse energy and multiplicity, elliptic flow); correlations and 
fluctuations: (fluctuations in charge and (pr), HBT measurements), hadron 
spectra: (low-pr single-hadron spectra and radial flow, particle ratios, reso- 
nances, anomalous p/n ratio at intermediate pr)', high-p T physics: (high-p T 
singles spectra, suppression phenomena in A + A, nonsuppression in d + A, 
high-pr two-particle correlations, nuclear suppression/enhancement in for- 
ward/backward directions), heavy flavor production: (charm, J/ip), and elec- 
tromagnetic probes: (direct photons). However, an important conclusion of 
this report is that systematic studies of these observables (vs. collision species 
and energy) are needed to extract unambiguous information on most of these 
features. 



1.3 Organization of this Document 

As a result, this paper concentrates on those aspects of the present data that 
address the broad features of energy density, thermalization, deconfinement 
and critical behavior. The focus in most cases will be on the data of the 
PHENIX experiment, but the data of the other RHIC experiments will be 
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cited to support and to extend the discussion 3 . The experimental tools that 
allow the systematic study of all phenomena as a function of the inferred im- 
pact parameter are presented in the context of hard-scattering phenomena. 
These methods and the associated data are then used to discuss the exper- 
imental evidence for the formation of a state of high-density matter. The 
measured abundances, spectra and flow patterns are used to analyze the de- 
gree of thermalization and collectivity in the produced matter. These results 
are then examined for evidence establishing the role of deconfined quarks and 
gluons in the produced system, along with the implications for its description 
as a quark-gluon plasma. A concluding section summarizes the findings and 
identifies key future measurements required to further refine our observations. 



2 ENERGY DENSITY AND E T , N CH 



A prerequisite for creating a quark-gluon plasma is producing a system with 
sufficiently large energy density. From both elementary estimates [22] and 
from extensive numerical studies in lattice QCD [11,10], the required density 
is known to be on the order of 1 GeV/fm 3 . Establishing that this energy 
density is created in RHIC collisions is a basic ingredient in establishing the 
creation of a QGP at RHIC. 

In this section we explore what can be deduced about the energy densities 
achieved in RHIC A + A collisions from measurements of the global transverse 
energy and multiplicity. In later sections these estimates will be compared to 
densities inferred from hydrodynamics-based models (Section 3) and from jet 
quenching evidence (Section 6). 

Specifically, we will address three different energy density estimates, and in- 
troduce two distinct time scales: (i) The peak general energy density that 
is achieved when the incoming nuclei overlap; (ii) The peak formed energy 
density involving created particles at proper time TForm', and (iii) The peak 
thermalized energy density present at proper time TTherm when local thermal 
equilibrium is first achieved (assuming that this occurs). The values and time 
scales for formed and thermalized energy densities are indicated schematically 
in Fig 4; detailed explanations follow in Sections 2.3 and 2.4. 

In this Section we will also review data on overall particle multiplicities, and 
through them distinguish between different models of the initial particle pro- 
duction. 



3 An underappreciated aspect of the RHIC program is the excellent agreement 
between the various experiments in almost all measured channels. 



17 




Fig. 4. Schematic drawing of the time and energy density scales derived through 
the Bjorken picture. 

2.1 General Energy Density 

The simplest definition of "energy density" is the total mass-energy within 
some region of space divided by the volume of that region, as seen at some 
instant of time in some Lorentz frame. However, this definition is not satisfac- 
tory since we can "trivially" raise any simple energy density by viewing the 
system in a different frame. For example, a static system with constant energy 
density po in its rest frame — say, a gold nucleus — will appear to have energy 
density 7 2 po when viewed in a frame boosted by Lorentz 7. Accordingly, we 
can only calculate a meaningful energy density (e) as mass-energy/volume for 
some region in the case when the total momentum in the region is zero. 

Now let us imagine a symmetric RHIC A + A collision at a moment when the 
two original nuclei are overlapping in space, as seen in the laboratory/center- 
of-mass frame. The total momentum in any overlap region is zero by symmetry, 
so we can calculate a meaningful — if short-lived — energy density for such a 
region. If each nucleus has energy density po i n its rest frame then the total 
energy density in the overlap region is just (e) = 2p 7 2 - If we take a nominal 
Po = 0.14 GeV/fm 3 for a nucleus at rest and 7 = 106 for a full-energy RHIC 
collision, then the result for the peak general energy density is (e) =3150 
GeV/fm 3 . This is a spectacularly, almost absurdly high number on the scale 
of ~1 GeV/fm 3 associated with the familiar transition described by lattice 
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QCD. 



This energy density is of course artificial, in that it would be temporarily 
present even in the case of no interactions between the two nuclei. It is in- 
structive to consider the (again artificial) case where the nucleons in the two 
nuclei have only elastic interactions. Then the time during which a high en- 
ergy density is present over any volume cannot last longer than t = 2R/ r y, 
where R is the rest-frame radius of the nucleus. With R =7 fm for Au this 
time is only 0.13 fm/c at RHIC, and after this time all energy densities will 
fall precipitously back to p if no secondary particles are created. The scale of 
this interval is so short that a scattering cannot even be said to have occurred 
within that volume unless its momentum transfer scale Q exceeds at least 1.5 
GeV/c, or more. Accordingly, we will turn our attention instead to energy 
densities involving only produced particles, as the potential source for a QCD 
transition. 

2.2 Formed Energy Density 

In any frame (not just the center-of-mass frame) where the two incoming nuclei 
have very high energies the region when/ where the nuclei overlap will be very 
thin in the longitudinal direction and very short in duration. In this limit, then, 
it is fair to describe all secondary produced particles as having been radiated 
out from a very thin "disk", and that they are all created at essentially the 
same time. These realizations lead directly to the picture described by Bjorken 
[74], whose original diagram is reproduced in Fig. 5 and whose derivation we 
retrace briefly here. 

Once the beam "pancakes" recede after their initial overlap, the region between 
them is occupied by secondaries at intermediate rapidities. We can calculate 
the local energy density of these created particles if we make one further 
assumption: that the secondaries can be considered "formed" at some proper 
time TForm after they are radiated out from the thin source disk. 

Our region of interest, in any frame, will be a slab perpendicular to the beam 
direction, with longitudinal thickness dz, with one face on the "source" plane 
in this frame, and transverse extent with area A covering the nuclear overlap 
region 4 . At time t = TForm this volume will contain all the (now- formed) par- 
ticles with longitudinal velocities < (3\\ < dz/ TForm (since we assume that the 
particles cannot scatter before they are formed!). We can then write this num- 
ber of particles as dN = (dz/T Fo rm)jj^, or equivalently dN = (dz/ ' TForm) ^r, 

4 The region described here corresponds to half the shaded region shown in Fig. 
5. Since /3y ~ for particles near the source location, this is an appropriate region 
over which we can calculate a meaningful energy density. 
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Fig. 5. Figure from Bjorken [74] illustrating the geometry of initially produced 
particles at a time t after the overlap of the incoming nuclei in some frame. The 
picture is valid in any frame in which the incoming nuclei have very high energies 
and so are highly Lorentz contracted. 



where y is longitudinal rapidity, since dy = d(3\\ at y — (3\\ = 0. If these particles 
have an average total energy (ijit) in this frame (E = m F for particles with 
no longitudinal velocity), then the total energy divided by the total volume of 
the slab at t — Tp orm is just 



dN(m T ) dN{r Form ) (m T ) 

{£{T Form )) = 



dz A dy T Form A 

1 dE T (r Form ) 

TFormA dy 
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where we have equated = (my)^ and emphasized that Eq. 5 is true for 
the transverse energy density present at time t = Tp orm . 

Equation 5 here is essentially identical 5 to Eq. 4 of Bjorken's result [74], 
and so is usually referred to as the Bjorken energy density e B j. It should be 
valid as a measure of peak energy density in created particles, on very general 
grounds and in all frames, as long as two conditions are satisfied: (1) A finite 
formation time TForm can meaningfully be defined for the created secondaries; 
and (2) The thickness/ "crossing time" of the source disk is small compared to 
TForm, that is, T Form ^> 2R/ r y. In particular, the validity of Eq. 5 is completely 
independent of the shape of the dE T {jF orm )l dy distribution to the extent that 
/3\\ is infinitesimally small in a comoving frame; a plateau in dE^/dy is not 
required. For present practical purposes, we will consider condition (2) above 
to be satisfied as long as r Form > 2R/ r y is true, corresponding to r Form >0.13 
fm/c for full-energy Au+Au collisions at RHIC. 

Bjorken's original motivation was to estimate, in advance of data, the energy 
densities that would be reached in high-energy A + A collisions, using knowl- 
edge of p(p) +p collisions to estimate (m T ) and dN/dy, and choosing Tp orm ~1 
fm/c without any particular justification other than as an order-of-magnitude 
estimate. With A + A collision data in hand, attempts have been made to use 
Eq. 5 to estimate the energy densities that are actually achieved in the colli- 
sions. Historically, e B j has been calculated using the final-state dE T /dy and 
simply inserting a nominal value of 1 fm/c for Tp orm . In addition, fixed target 
experiments have been using dEp/dr] as an estimate for dEp/dy, which is a 
good approximation for these experiments; at RHIC a correction is made for 
the Jacobian dy/dr] which is important for a collider geometry. These "nominal 
Bjorken energy density" estimates, which we term e I ^° mmal ^ range for central 
event samples from about 1.5 GeV/fm 3 in Au+Au collisions at AGS energies 
I 75 ] = 5 GeV )> t0 about 2 - 9 GeV/fm 3 in Pb+Pb collisions at SPS en- 

ergies [76,38] (y/s~^ =17 GeV; and see also [59]) to about 5.4 GeV/fm 3 in 
Au+Au collisions at full RHIC energy [59] ( v /s7^7 =200 GeV). 

It has often been noted that all of these values are similar to, or higher than, 
the 1 GeV/fm 3 scale required for the QCD transition. However, we cannot 
take these £^° mmal estimates seriously as produced energy densities without 
some justification for the value of 1 fm/c taken for TForm- An indication of 
potential problems with this choice arises immediately when considering AGS 
Au+Au and SPS Pb+Pb collisions, where the center-of-mass "crossing times" 
2i?/7 are 5.3 fm/c and 1.6 fm/c respectively, which implies that this choice for 
TForm =1 fm/c actually violates the validity condition TForm > 2R/l we set 
for the use of Eq. 5. So we will deprecate the use of £^? minal ag an q Uan titative 
estimate of actual produced energy density, and instead treat it only as a 



A (well-known) factor of 2 error appears in the original. 
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compact way of comparing dEr/dr] measurements across different systems, 
centralities and beam energies. 



2.3 Realistic r Form and Ebj estimates 



Can we justify a better estimate for T Form ? We might say, on general quantum 
mechanical grounds, that in a frame where its motion is entirely transverse 
a particle of energy mx can be considered to have "formed" after a time 
t = h/rriT since its creation in that frame. To estimate the average transverse 
mass, we can use the final-state dE F /dr] to estimate dE F {r Form )j dy and, cor- 
respondingly, use the final-state dN/drj as an estimate for dN{r Form ) / dy to 
obtain 

dE T (T Form ) / dy dE T /dr) 

\ m T) = —tttt rh— ^ —rrrn— (Final state). (6) 

dN(T Form )/dy dN/dr) 



PHENIX has measured the ratio of final-state transverse-energy density to 
charged-particle density, each per unit pseudorapidity, and the results are 
shown in Fig. 6. For a wide range of centralities the ratio is remarkably con- 
stant at about 0.85 GeV for full-energy central Au+Au collisions and shows 
very little change with beam energy, decreasing to only 0.7 GeV when ^/s NN 
is decreased by an order of magnitude down to 19.6 GeV. 

If we approximate dN Ch /dr/ = (2/3)dN/drj in the final state then Eq. 6 would 
imply (m T ) ~ 0.57 GeV and corresponding r Form ~ 0.35 fm/c, a value shorter 
than the "nominal" 1 fm/c but still long enough to satisfy our validity condi- 
tion T Form > 2/2/7 & t RHIC. Inserting this value into Eq. 5, along with the 
highest dEr/dr] = 600 GeV for 0-5% central events as measured by PHENIX 
[59], yields a value of (e) = 15 GeV/fm 3 for the energy density in initially 
produced, mid-rapidity particles in a central RHIC Au+Au collision, that is, 
roughly 100 times the mass-energy density of cold nuclear matter. 

It is important to note that this large value of the energy density as ob- 
tained from Eq. 5 represents a conservative lower limit on the actual (s(T Form )) 
achieved in RHIC collisions. This follows from two observations: (1) The final- 
state measured dEx/drj is a solid lower limit on the dET(r Form ) / dy present 
at formation time; and (2) The final-state ratio (dE T / drj) / (dN / drj) is a good 
lower limit on (mr) at formation time, and so yields a good upper limit on 
TForm- We justify these statements as follows: 

Several mechanisms are known that will decrease dE F /dy as the collision sys- 
tem evolves after the initial particle formation, while no mechanism is known 
that can cause it to increase (for y — 0, at least). Therefore, its final-state 
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Fig. 6. The ratio of transverse energy density in pseudorapidity to charged particle 
density in pseudorapidity, at mid-rapidity; shown as a function of centrality, rep- 
resented by the number of nucleons participating in the collision, N par t, for three 
different RHIC beam energies [59]. 

value should be a solid lower limit on its value at any earlier time. A partial 
list of the mechanisms through which dE^jdy will decrease after t = Tp orm 
includes: (i) The initially formed secondaries in any local transverse "slab" 
will, in a comoving frame, have all their energy in transverse motion and none 
in longitudinal motion; if they start to collide and thermalize, at least some 
of their E? will be converted to longitudinal modes in the local frame; (ii) 
Should rough local thermal equilibrium be obtained while the system's ex- 
pansion is still primarily longitudinal, then each local fluid element will lose 
internal energy through pdV work and so its Et will decrease; (iii) If there are 
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pressure gradients during a longitudinal hydrodynamic expansion then some 
fluid elements may be accelerated to higher or lower rapidities; these effects 
are complicated to predict, but we can state generally that they will always 
tend to decrease dE T /dy where it has its maximum, namely at y — 0. Given 
that we have strong evidence that thermalization and hydrodynamical evolu- 
tion do occur in RHIC collisions (Section 3), it is likely that all these effects 
are present to some degree, and so we should suspect that final-state dEx/dr) 
is substantially lower than dE T {jF orm )l dy at mid-rapidity. 

Turning to our estimate of TForm, the assumption that TForm = h/ (mr) cannot 
be taken as exact, even if the produced particles' m^'s are all identical, since 
"formed" is not an exact concept. However, if we accept the basic validity of 
this uncertainty principle argument, then we can see that the approximation 
in Eq. 6 provides a lower limit on (m T ). First, the numerator dE T /dr] is a 
lower limit on dE T (TF orm )/dy, as above. Second, the argument is often made 
on grounds of entropy conservation that the local number density of particles 
can never decrease [77], which would make the final-state denominator in Eq. 
6 an upper limit on its early-time value. 

With these limits in mind, then, it is not unreasonable for us to claim that the 
peak energy density of created particles reached in central Au+Au collision 
at RHIC is at least 15 GeV/fm 3 , and in all likelihood is significantly higher. 

2.4 Thermalized Energy Density 

We have arrived at a reasonably solid, lower-limit estimate for the energy den- 
sity in produced particles in a RHIC Au+Au collision, and it is more than 
enough to drive a QCD transition. But the situation at t — TForm pictured in 
Fig. 5 looks nothing like local thermal equilibrium. It is an important ques- 
tion, then, to ask: if and when the system evolves to a state of local thermal 
equilibrium, is the energy density still sufficient to drive the transition to a 
QGP? 

To answer this we begin by looking at the state of the system at t — TForm 
and immediately afterward. At the time they are formed the particles have 
sorted themselves out automatically, with all the particles on a "sheet" at a 
longitudinal position z having the same longitudinal velocity f3\\ = z/t; and 
so in the rest frame of a sheet all the sheet's particles have only transverse 
motion. If the particles continue free-streaming and never reinteract then the 
energy density will continue to fall as e ~ 1/t and the Bjorken formula in Eq. 

5 will be valid, with t in place of TForm, as long as the expansion is primarily 
longitudinal 6 . 

6 For long times t > R transverse expansion will become significant and the energy 
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For thermalization to occur the particles will have to start interacting and/or 
radiating. Once this happens the particles which were originally together on 
one "sheet" will start to spread in longitudinal velocity, though on short time 
scales we would expect their group average longitudinal velocity to remain the 
same. If the thermalization process is fast enough, then, we would expect that 
at time t = TTherm these groups will have formed locally equilibrated fluid 
elements, with a velocity profile following f3P uld = z/t. The energy density at 
this time will be reduced from the energy density at formation time e(TForm) 
by a factor T Form /T Therm ; i.e. the Ssj of Eq. 5 but with r Therm in place of 
r Form . This evolution is illustrated in Fig. 4. 

Once local equilibration is achieved we would then expect the system to evolve 
hydrodynamically, and the behavior of e(t) will depend on the details of the 
local equations of state (EOS). Without knowing those details, though, we can 
say that in the limit of low pressure, p/e ~ 0, the energy density will continue 
to evolve (during longitudinal expansion) as e ~ 1/t, while in the limit of 
high pressure, p/e ~ 1/3, the energy density will decrease somewhat more 
quickly, e ~ l/£ 4//3 , within a fluid element. This range of possible behaviors 
for t > TTherm is indicated schematically in Fig. 4. 

A direct theoretical determination of TTherm would require a detailed descrip- 
tion of both the parton-parton interactions and the resulting evolution of the 
system density. However, other lines of reasoning may provide information 
on TTherm- For example, it has been argued [78] that the strong elliptic flow 
in RHIC collisions can be taken as evidence for fast thermalization (see Sec- 
tion 3.3). In a hydrodynamic picture the source of elliptic flow is the spatial 
anisotropy of the energy density in the transverse plane at the time hydrody- 
namics becomes valid. If local equilibration and the onset of hydrodynamics 
is delayed because interactions between the initially produced particles are 
weak at first, then the spatial anisotropy which could give rise to elliptic flow 
will be reduced (see Fig. 14). This, in effect, limits how high TTherm can be if 
hydrodynamics is the mechanism for generating elliptic flow. 

We can see from Table 1 in Section 3.5 that hydrodynamical models typically 
require quite short thermalization times, in the range of 0.6-1.0 fm/c, in order 
to reproduce the magnitude of elliptic flow which is observed at RHIC. If we 
take this range as typical of what hydrodynamics would imply for TTherm, then 
we can calculate the corresponding "typical" implied energy densities at ther- 
malization time as in range of 5.4 GeV/fm 3 to 9.0 GeV/fm 3 . These densities 
are well above that required to drive the QCD transition, so the combination 
of our transverse energy measurements and the fast thermalization times from 
hydrodynamics can be taken, to some degree, as evidence that conditions to 
create the equilibrated upper phase of QCD matter are achieved at RHIC. 



density will decrease as e ~ 1/t 3 . 
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2.5 What Are the Initial Quanta? 



With our extensive use of the picture in Fig. 5 it is only natural to ask, "What 
are these initially produced particles?" that Bjorken referred to, nonspecifi- 
cally, as "quanta" . What models do we have for initial production, and what 
can we say about them using our data on Et and multiplicity? 

The simplest assumption is that the initially produced particles in a RHIC 
collision are scattered partons at mid- to low-p^, traditionally known as "mini- 
jets" . For a long period in advance of RHIC data, it was widely expected that 
mini-jets would be the dominant channel for Et and particle production, and 
this led to two further, general expectations: first, that multiplicity and Et 
per interacting nucleon would go up sharply at collider energies, as compared 
to fixed-target energies, since jet and mini-jet cross sections are increasing 
quickly with (see Fig. 7); and, secondly, that E T and multiplicity per 
participating nucleon would increase steeply in more central events, since the 
rate of hard pQCD scatterings goes up faster with centrality than does the 
number of interacting nucleons. 

It was therefore quite surprising when the first RHIC data [80,37,38] showed 
lower multiplicities than had been predicted from mini-jet models, and only a 
modest increase in Et and multiplicity per participant as functions of central- 
ity. Compared to the sharp rise, shown in Fig. 7, predicted by straightforward 
factorized pQCD, it was clear that some mechanism must be acting at RHIC 
energies to restrict, or regulate, particle production [81,79]. 

pQCD-based models have parameters regulating the momentum scales; these 
include a lower-momentum cutoff, and the factorization and fragmentation 
scales. Figure 8 shows that the pQCD-based HIJING model, circa 2002, was 
able to reproduce 130 GeV and 200 GeV dN c h/dr] reasonably well. However, 
in that model jet production via hard scattering is an important mechanism 
for particle production, and the combination of the ^/s dependence of hard- 
scattering cross sections with the growth of the nuclear overlap with centrality 
causes the model to predict an increase in the ratio between the two data sets 
with centrality. The observed ratio is, instead, quite constant. Thus the authors 
found it necessary to introduce a centrality-dependent shadowing to regulate 
the jet growth [79]. 

An alternative to models which use collinearly factorized pQCD is found in 
the "color glass condensate" picture, in which the gluon population of low- 
x, low-pT states in the initial nuclear wave function is limited by transverse 
overlap and fusion of these low-pT gluons. The phase-space density saturates 
because of the competition between extra gluon radiation from higher-x gluons 
and non-linear fusion of the gluons at high density. Au+Au collisions are then 
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Fig. 7. Figure from Li and Wang [79] showing trends in final-state charged mul- 
tiplicity per participant pair vs. (nucleon-nucleon) beam energy. (PHENIX data 
points[59] have been added.) The curves are the result of their two-component 
"hard/soft" model, which reproduces well the multiplicities from elementary p{p)+p 
collisions at RHIC energies. The same model extended to nuclear collisions with no 
regulating mechanism on hard processes (the "No Shadowing" line) over-predicts 
the multiplicities in central RHIC collisions, while the data can be matched if sub- 
stantial nuclear shadowing of gluons is invoked (shaded band) . 

collisions of two sheets of colored glass, with the produced quarks and gluons 
materializing at a time given by the inverse of the saturation momentum, 
r = 1/Qs- Saturation of gluons with momenta below Q s provides a regulating 
mechanism that limits the rise in gluon — and later, hadron — multiplicity with 
centrality and beam energy. Models featuring this initial-state gluon saturation 
agree well with essentially all RHIC data on the multiplicity density, which is 
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Fig. 8. Multiplicity per participant nucleoli pair, as a function of centrality, for 
y/s^: =130 GeV and 200 GeV Au+Au collisions as measured in PHENIX [59]; com- 
pared to theoretical predictions available in 2002. "HIJING" is a pQCD-based model 
[82], while "KLN" features gluon saturation in the initial state [83,84]; "EKRT" as- 
sumes saturation in the final state [85,86]. 

dominated by low-momentum particles [83,84]. This is seen, for instance, in 
Fig. 8. 

In this picture, the total gluon multiplicity is proportional to l/a s • Q 2 S1 which 
limits the number of low-momentum charged particles produced. Q s evolves 
slowly with collision centrality and beam energy. For central Au+Au collisions, 
it has been estimated that the typical m T scale of the gluons "liberated" from 
the colored glass is about 1 GeV per particle [77] , which is above the lower limit 
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of 0.53 GeV per particle that we set above using the PHENIX data. Though 
there are fewer predictions of E? than total charged-particle production from 
gluon-saturation models, the existing models are broadly consistent with data 
at RHIC. Consequently, gluon saturation is considered to be a promising can- 
didate for describing the initial state of RHIC collisions. 

2. 6 Conclusions 

Using reasoning similar to that of Bjorken [74], combined with some simple 
formation-time arguments, we can draw the following conclusions from the 
PHENIX data on transverse energy production and overall particle multiplic- 
ity: 

• The peak energy density in created secondary particles is at least 15 GeV/fm 3 , 
and this is most likely an underestimate. This is well in excess of the ~1 
GeV/fm 3 required, according to lattice QCD predictions, to drive a QCD 
transition to QGP. 

• We note that hydrodynamical calculations which reproduce the magnitude 
of elliptic flow observed at RHIC require local thermalization to occur very 
quickly, typically by 1 fm/c or earlier (see Section 3.5). If the system does 
reach local equilibrium on this time scale then the energy density of the first 
thermalized state would be in excess of 5 GeV/fm 3 , well above the amount 
required to create the QGP. 

• Pre-RHIC expectations that E T and charged particle production would be 
dominated by factorized pQCD processes were contradicted by data, which 
showed only very modest increases with centrality and beam energy. A new 
class of models featuring initial-state gluon saturation compares well with 
RHIC multiplicity and Et data, and are also consistent with our Bjorken- 
style arguments for estimating energy densities at early times. 



3 THERMALIZATION 

A key question is whether the matter formed at RHIC is thermalized, and if 
so when in the collision was equilibration achieved. If thermalization is estab- 
lished early then evidence for strong transverse expansion can be potentially 
related to the equation of state of the dense matter produced at RHIC. To 
explore these issues we review several experimental observables from integral 
quantities (numbers of particles produced and in what ratios), to differen- 
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tial distributions (measured px and azimuthal distributions), to two-particle 
(HBT) correlations. 

3. 1 Chemical Equilibrium 

For many years it has been known that the abundances of different hadron 
species in e + + e~ and p + p reactions can be reproduced by statistical models 
[87,88]. This success is often attributed to hadronization statistically filling the 
available phase space. At RHIC there is also the possibility that the strong 
scattering deduced from the measurements of elliptic flow (section 3.3) may 
prove sufficient to establish chemical equilibrium. 

The production of strange particles provides a means to check whether chem- 
ical equilibrium is achieved. For e + + e~ and p + p reactions strange particle 
production is suppressed due to the small size of the system. This canonical 
suppression is largely removed for central heavy-ion collisions. If the mea- 
sured strangeness yields are still lower than full equilibrium predictions, then 
the partial equilibrium can be quantified by a multiplicative factor of 7 S for 
each strange quark in a hadron, where 7 S = 1 for complete equilibration and 
7 S < 1 for partial equilibration. 

Figure 9 shows the centrality dependence of K/n and p/n ratios in Au+Au 
collisions at ^/sntv — 200 GeV [54]. Both K + /tt + and K~ /n~ increase rapidly 
for peripheral collisions, and then saturate or rise slowly from mid-central to 
the most central collisions. The ratios p/n + and p/n" also increase from pe- 
ripheral collisions but appear flatter than the Kj n ratios. Canonical statistical 
models [89] predict an increase in these ratios with centrality, as the larger 
system-size effectively places less of a constraint on conserved quantities. In 
addition the chemical parameters, T c hem and /ib, can also vary with centrality 
[90,91]. 

Focusing on the ratios from central collisions at y^J/vjv = 200 GeV, the data 
are compared to the thermal model data analysis of Kaneta and Xu [90] in 
Fig. 10. The extracted thermal parameters from this fit are T chem = 157 ± 3 
MeV, \i B = 23 ± 3 MeV, and 7 S = 1.03 ± 0.04. A large 7 S is also found by 
STAR [92] who extract 7 S = 0.96 ±0.06, while Cleymans et al. [91] extract 7, 
that increases from 7 S ~ 0.85 in peripheral collisions to 7 S ~ 0.95 for central 
collisions at RHIC. Similar fits to the central RHIC data are obtained by 
Braun-Munzinger et al. [93] who assume complete chemical equilibration, i.e. 

Is = 1- 

We note that there are differences in the temperature parameter extracted by 
the different authors. Kaneta and Xu [90] extract T chem = 157 ±3 MeV which 
is lower than that extracted by both Braun-Munzinger et al. [89] of T c h em = 
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Fig. 9. Centrality dependence of particle ratios for (a) K + /ir + , (b) K /it , (c) 
p/ir + , and (d) p/tt~ in Au+Au collisions at ^/snn = 200 GeV [54]. 



177 ± 7 MeV and Cleymans al. [91] of T chem = 165 ± 7 MeV. However, both 
Braun-Munzinger et al. [89] and Magestro [94] discuss the sensitivity of the 
extracted temperature to corrections from feed-down from decays. Cleymans 
et al. [91] estimate that over 70% of ir + in the thermal model fits come from 
the decay of resonances. 

At lower beam energies there is controversy over whether strangeness is in 
full chemical equilibrium. Becattini et al. [95] use data that is integrated over 
the full rapidity and find that strangeness is in partial equilibrium, i. e. at the 
AGS 7, = 0.65 ±0.07 and at the SPS 7, = 0.84 ±0.03. Braun-Munzinger et al. 
[89] instead use ratios measured at mid-rapidity which typically have larger 
strange/no- nstrange values, and, hence, they obtain acceptable fits with 7 S = 1 
at both AGS and SPS energies. At RHIC energies thermal model comparisons 
all use mid-rapidity data; a choice that is motivated in part by the separation 
between fragmentation regions and central particle production. 
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Fig. 10. Comparison of PHENIX (triangles), STAR (stars), BRAHMS (circles), and 
PHOBOS (crosses) particle ratios from central Au+Au collisions at ^/snn = 200 
GeV at mid-rapidity. The thermal model descriptions from Kaneta [90] are also 
shown as lines. See Kaneta [90] for the experimental references. 

In contrast to the controversies at lower beam energies, the observation that 
strangeness is equilibrated is common to all thermal calculations that repro- 
duce RHIC data. This is consistent with chemical equilibrium being obtained 
before hadronization, though does not prove that this is the case. An alter- 
native explanation is that scattering in the hadronic phase could increase 7 S 
to 1, though small interaction cross sections imply that it may be difficult to 
equilibrate the multi- strange baryons before the hadrons freezeout. 



3. 2 Spectra 



Hadron spectra reflect conditions late in the reaction, as well as the integrated 
effects of expansion from the beginning of the collision. Figure 11 shows the pr 
distributions for pions, kaons, protons, and anti-protons in both central (top 
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panel) and peripheral collisions (bottom panel) [54]. The pion spectra have a 
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Fig. 11. Transverse momentum distributions for pions, kaons, protons, and 
anti-protons in Au+Au collisions at y^sjvTv = 200 GeV [54]. 

concave shape at low pr where many of the pions may come from the decay of 
resonances: A, p etc. The kaon spectra are approximately exponential over the 
full measured p^ range, whereas the proton spectra flatten at low px for the 
most central collisions. A striking feature is that the proton and anti-proton 
spectra in central collisions become comparable in yield to the pion spectra 
above 2 GeV/c. This is more fully discussed in Section 7. 

One way to characterize the change in spectra as a function of centrality 
is to calculate (p T ) for each spectrum [54] as shown in Fig. 12. The (p T ) 
increases for all particles as a function of centrality with the largest change 
occurring in peripheral collisions (N part < 100). Across the different particles 
the increase is largest for protons and anti-protons. This is consistent with 
a collective expansion velocity that increases with centrality to produce the 
largest increase in (p T ) for the heaviest particles. 

The pion, kaon, and proton spectra can all be fit using an ansatz of a thermal, 
expanding source [96,48] to extract the collective transverse expansion velocity 
(Pt) as well as the temperature at freezeout, Tj D . Figure 13 shows (Pt) ~ 0.45 
at AGS energies [97,98], which increases to (f3 T ) ~ 0.5 at the SPS [99,100,101] 
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Fig. 12. Mean transverse momentum as a function of AT par . t for pions, kaons, protons 
and anti-protons in Au+Au collisions at t/snn = 200 GeV [54]. The systematic 
errors from extrapolation, which are scaled by a factor of two for clarity, are shown in 
the bottom for protons and anti-protons (dashed-dotted lines), kaons (dotted lines), 
and pions (dashed lines). The shaded bars to the right represent the systematic error. 

and RHIC [48,102]. All the above fits use similar model assumptions of a linear 
velocity profile and a Woods-Saxon density profile. That the spectra at these 
beam energies can be reproduced by a thermal source is necessary, but not 
sufficient, evidence for thermal equilibrium at each of these energies. However, 
it is difficult to draw strong conclusions from the increase in (/3r) as a function 
of beam energy since the parameters (fir) and Tf Q are strongly anticorrelated 
and their values depend on fit ranges and treatment of decays. 



3.3 Elliptic Flow 



At the beginning of a heavy ion collision, the spatial distribution of the col- 
liding matter resembles an ellipsoid due to the incomplete overlap of the two 
colliding nuclei. Any strong scattering in this early stage converts the spatial 
anisotropy to a momentum anisotropy which is observable as an elliptic flow of 
the emitted hadrons. Elliptic flow is a self-limiting phenomenon, which is read- 
ily understood in the thermodynamic limit. If strong scattering is sufficient 
to establish local thermal equilibrium, then the pressure gradient is largest in 
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Fig. 13. Beam-energy dependence of the extracted mean transverse expansion ve- 
locity as a function of beam energy from simultaneous fits to spectra of different 
mass [97,98,99,100,101,48,102]. 

the shortest direction of the ellipsoid. This gradient produces higher momenta 
in that direction, quickly reducing the spatial asymmetry. 

The absence of any strong scattering in the early stage of the colliison would 
reduce the amount of elliptic flow that could be created. If the initially pro- 
duced particles are allowed to initially free stream and reach local equilibrium 
only after some time delay, then the spatial anisotropy at the start of hy- 
drodynamic evolution will be reduced; the longer the delay, the greater the 
reduction. Following the prescription of Kolb et al. [78] , we plot in Fig. 14 the 
eccentricity after a time delay At compared to its value at formation time, as 
a function of Au+Au collision centrality. The eccentricity (e) of the reaction 
zone is 



(y 2 ) - (x 2 ) 
(y 2 ) + (x 2 )' 



(7) 



The eccentricity can be analytically calculated once the density profile of the 
nuclei is chosen (typically a Woods-Saxon shape). It can also be calculated 
using Monte Carlo techniques, where the positions of those nucleons that 
participate in the reaction are used to calculate the averages in Eq. 7. From 
Fig. 14 we can see that for time delays of 2 fm/c or greater the magnitude 
of the eccentricity is significantly reduced, and its shape vs. centrality is also 
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altered. 



If locally equilibrated hydrodynamics is taken as the mechanism for generating 
elliptic flow, then the observation of any substantial amount of elliptic flow 
can be taken as evidence that local thermal equilibrium is achieved on a time 
scale before the spatial anisotropy would be completely erased. The order of 
this time-scale would be t ~ R/c, where R is the nuclear radius. However, 
the hydrodynamical calculations we will examine here (see Sec. 3.5 and Table 
1) all require quite short thermalization times, from 0.6-1.0 fm/c, in order to 
reproduce the magnitude of elliptic flow observed at RHIC. 
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Fig. 14. The ratio of the eccentricity after a time delay At compared to its value at 
formation time, as a function of Au+Au collision centrality. The calculations follow 
the prescription of [78] where the produced particles are allowed to free-stream at 
first and reach local equilibrium only after some time delay. 



The azimuthal anisotropy of the spectra can be characterized in terms of 
Fourier coefficients, which at RHIC are dominated by the elliptic flow, the 
second Fourier coefficient, i^Cpt), where 

d 2 N 

N (l + 2v 2 ( PT ) cos(20)) . (8) 



Both the first Fourier coefficient, vi, and higher order coefficients have been 
neglected in the above expression. 
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The most direct evidence that i>2 is related to spatial asymmetries present 
early in the reaction is that v 2 at low approximately scales with the initial 
eccentricity (e) of the reaction zone. The measured values of v-i normalized by 
e are shown in Fig. 15 vs. centrality for two different p T ranges [46]. At low 
momentum t> 2 /e is independent of centrality to within 20%. This scaling is 
increasingly broken at higher pt- 
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Fig. 15. A2 = V2/E vs. centrality for Au+Au collisions at y / ijv r /v = 130 GeV [46]. 
The data points come from two different types of two-particle correlations: "fixed" 
Pt correlations when both particles are at the same pr (points are labeled as "F"), 
and "assorted" pr correlations when the two particles have different pt (points 
are labeled as "A"). In this case the labeled pt range is for the higher-momentum 
particle of the pair. 



The measured values of the integrated v 2 at RHIC are larger than those at 
lower energies, but this is in part due to the fact that ^(pt) increases with 
Pt and (pt) increases as a function of beam energy. To remove this effect we 
will concentrate on the differential flow, i.e. the shape of ^(pt)- 

To make a uniform comparison between different colliding nuclei (Pb+Pb at 
SPS and Au+Au at RHIC) as well as different impact parameter selections 
from the different experiments, we normalize v 2 by the eccentricity, e, as shown 
in Fig. 16. The values of e have been calculated via a Glauber Monte Carlo 
using Woods-Saxon density distributions for the Au and Pb nuclei. The aver- 



37 



w 

> 



1 

0.8 
0.6 
0.4 
0.2 




• 7i", K", PHENIX, 200 GeV, 20-40% 

■ n, STAR, 1 30 GeV, 1 1 -46% 

▲ h*, PHENIX, 62 GeV, 20-40% 

O 7i, NA49, 17 GeV, 13-34% 

O 7C, CERES, 17 GeV, 13-26% 

□ 7t, NA49, 8 GeV, 1 3-34% 

• 4 



□ 



t 



o 




8 



0.5 



1.5 



T 

ft 
1 



2 2.5 

p T (GeV/cj 



Fig. 16. V2{pt)/z vs. j>t for mid-central collisions at RHIC (filled symbols) and SPS 
(open symbols). Dividing by eccentricity removes to first order the effect of different 
centrality selections across the experiments [50,103,70,104,105]. 



ages in Eq. 7 are over the participating nucleons, hence e is calculated at the 
start of the collision. The pion data in Fig. 16 show that ^(pt)/^ increases 
approximately linearly with px for low Pt- The rate of increase of v-ije as a 
function of pt is larger at RHIC [50,103] than at SPS [104,105] as can most 
easily be seen by calculating the slope of v^je below pt =1 GeV/c (Fig. 17). 
The slope (dv 2 /dpT)/£ increases from SPS to RHIC by approximately 50%. 
Hydrodynamical calculations [106] shown in Fig. 17 reproduce the data both 
at RHIC and at CERN SPS within one standard deviation. More extensive 
comparisons with hydro calculations will be discussed in section 3.5, while the 
behavior of v 2 at higher px, which follows a scaling with respect to the number 
of quarks, is discussed in Section 7. 

Further insight into the expansion dynamics can be obtained from the mass 
dependence of ^(pr) shown in Fig. 18 for pions, kaons and protons [50] along 
with a comparison with an early hydrodynamic model calculation [107]. The 
V2{pt) for pions is larger than for kaons and protons at low px, and this mass 
ordering has been explained as resulting from radial expansion[107] that pro- 
duces a larger distortion of the elliplic flow induced velocity profile for larger 
hadron masses. However, as will be discussed in Section 3.5, this calculation 
fails to reproduce the proton spectra, and attempts to remedy this failure 
lead to calculations that no longer reproduce the measured v 2 for pions and 
protons. 
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Fig. 17. The slope of the scaled elliptic flow, (dv2/dpr)/s, for mid-central collisions 
at RHIC (filled symbols) and the SPS (open symbols). The slope is calculated from 
the data in Fig. 16 for the data px < 1 GeV/c. The solid error bars represent the 
total systematic error including the systematic error on v 2 and e [50,103,70,104]. 

3.4 HBT 



Bose-Einstein correlations between identical particles provide a measure of 
the space-time extent of the source at the end of the reaction. Because the 
extracted source parameters as measured by the HBT technique are driven 
by space-time correlations, HBT results are sensitive to expansion dynamics 
integrated throughout the collision. HBT measurements were originally moti- 
vated by theoretical predictions of a large source size and/or a long duration 
of particle emission [108,109,110] — which would result from the presence of a 
long-lived mixture of phases in the matter as it undergoes a first-order phase 
transition from a quark-gluon plasma back to the hadronic phase. 

In HBT analyses, multidimensional Gaussian fits are made to the normal- 
ized relative momentum distributions yielding fit parameters, -Ri on g, -Rside, -Rout 
[111], also referred to as HBT radii, where 

C 2 = 1 + Aexp(-i?g ide g s 2 ide - R 2 out q 2 ovlt - R\ ong qf ong )- (9) 



The coordinate system is chosen so that the longitudinal direction is parallel to 
the beam axis, the out direction is in the direction of the pair's total transverse 
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Fig. 18. V2(pt) for pions, kaons and protons produced in minimum-bias collisions 
at RHIC [50] compared to hydro calculations from Huovinen et al. [107]. 

momentum, and the side direction is in the transverse plane perpendicular to 
the out axis. For dynamic (i.e. expanding) sources, the HBT radii depend on 
the mean transverse momentum of the particle pairs, kx = \P\t + P2t|/2, 
and correspond to lengths of homogeneity: regions of the source which emit 
particles of similar momentum [112]. Measuring the kx dependence of HBT 
radii provides essential constraints on dynamical models [113]. In particular, 
the ratio -R ut/-R S ide is predicted to be larger than unity for sources which emit 
particles over a long time. 

The measured kr dependence of all radii [56] and the ratio R out /Rside are 
shown in Fig. 19, along with STAR results [114]. The data from PHENIX and 
STAR are in excellent agreement. Both sets of data have been corrected for 
Coulomb repulsion between the detected particles. 

The measured radii all decrease with increasing kx as expected for a rapidly 
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Fig. 19. The kx dependence of the Bertsch-Pratt parameters for tt + tt + (blue square) 
and tt~tt~ (red circle) for 0—30% centrality with statistical error bars and systematic 
error bands. Results from PHENIX [56] and STAR [114] are overlaid. 

expanding source. The ratio R ou t/Rsidc was measured to be 1 within errors, 
with a slight systematic decrease for increasing kr- As is discussed in the 
next section, these data have excluded the validity of a large majority of 
hydrodynamical models developed to describe Au+Au collisions at RHIC, 
indicating that in their present form these models do not describe well the 
space-time evolution of the Au+Au collisions. 



3.5 Hydrodynamic Model Comparisons 



Many of the experimental features in the spectra and elliptic flow are con- 
sistent with equilibrium being established early in the collision with large 
pressure gradients that drive a strong expansion. Moving from a statement 
of "consistency" to a statement that equilibrium has been "established" is 
difficult. Some progress can be made by comparing the data to hydrodynamic 
models that assume full equilibrium early in the collision. 

A variety of hydrodynamic models have been published. Our approach is to 
confront these models with the following broad set of data; ^(pr), spectra, 
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and HBT. In this paper we will not compare the data with hydro-inspired 
parameterized fits, e.g. blast-wave [115] or Buda-Lund [116] models, but will 
restrict ourselves to dynamical hydro models. 

In Figs. 20 and 21, hydro calculations that include a phase transition from 
the QGP phase to a hadronic phase are shown with solid lines, while hydro 
calculations that do not include a pure QGP phase at any stage in the dy- 
namics are drawn with dashed lines. The four calculations that include a QGP 
phase all assume an ideal gas EOS for the QGP phase, a resonance gas for the 
hadronic phase and connect the two using a first-order phase transition and 
a Maxwell construction. These calculations use latent heats that range from 
0.8 GeV/fm 3 (Teaney et al. [106]) to 1.15 GeV/fm 3 (Huovinen et al. [107] and 
Kolb et al. [117]), to 1.7 GeV/fm 3 (Hirano et al. [118,119]). For comparison 
the bag model of the nucleon with external bag pressure B = (230 MeV) 4 
and a T crit = 164 MeV produces a latent heat of 1.15 GeV/fm 3 [120]. The 
calculations that do not include a QGP phase (dashed lines) either include a 
hadron phase and a phase mixture by forcing the latent heat of the transition 
to infinity [106], or use an hadronic resonance gas equation of state, i.e. no 
mixed or QGP phases [107]. 

The calculations also differ in how they solve the hydro equations and how they 
treat the final hadronic phase. The work of Hirano, Tsuda, and Nara cited here 
are the only calculations in this paper that solve the hydro equations in 3D 
[118,119]. For the final hadronic stage Teaney [106] uses a hybrid model that 
couples the hadronic phase to RQMD to allow hadrons to freezeout according 
to their cross section, i. e. for chemical equilibrium to be broken in the hadronic 
phase. Hirano [118] and Kolb [120] both allow for partial chemical equilibrium 
by chemically freezing out earlier than the kinetic freezeout. This has been 
done in order to reproduce the large proton yield measured at RHIC (see later 
in this section). In contrast, Huovinen [107] maintains full chemical equilibrium 
throughout the hadronic phase. 

Figure 20 compares these calculations to the measured minimun-bias proton 
and pion Vi(^r)je. Minimim-bias results were chosen in order to have the 
broadest set of data and model calculations for comparison. The four calcula- 
tions that include a phase transition from the QGP phase to a hadronic phase 
(solid lines) reproduce the low-p^ proton data better than the two hydro cal- 
culations that do not have a QGP phase at any stage in the dynamics (dashed 
lines). The presence of the first-order QGP phase transition softens the EOS 
which reduces the elliptic flow. At higher p T there is considerable variation 
between the models. Part of this is due to how the final hadronic stage is 
modeled. For example, Kolb's (solid light-blue line) and Hirano's (solid dark- 
blue line) calculations allow for partial chemical equilibrium in the final stage 
compared to Huovinen (solid green line) which chemically freezes out late in 
the collision. The difference is observable above p T ~ 1 GeV/c. 
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Fig. 20. Top two panels: On the left, proton -V2(pt) vs. Pt for minimum-bias 
collisions at RHIC [50,103] are compared with hydro calculations [106,118,120,107], 
and on the right is the same comparison for pions. Bottom two panels: On the left, 
proton spectra for 0-5% collisions at RHIC [54] are compared with the same hydro 
calculations and on the the right is the same comparison for pions. 

The same hydro models are compared to the pion v 2{pt)/z measurements from 
STAR and PHENIX in Fig. 20. The Kolb (solid light-blue line) and Hirano 
(solid dark-blue line) calculations fail completely by predicting too strong a t>2. 
These two models have very similar partial chemical equilibrium assumptions 
in the late hadronic stage. It is worth noting that the Kolb calculation is the 
same as the Huovinen (solid green line) calculation with the exception of the 
final hadronic stage. 

All the above models have assumed ideal hydrodynamics, i. e. with no viscosity 
and zero mean free path. Non-zero viscosity in the QGP reduces t>2 [121,122] 
and since the early hydro calculations from Teaney and Huovinen reproduced 
the magnitude of the pion t> 2 data, it is often stated that viscosity of the 
matter at RHIC must be small [24] . However recent calculations from Hirano 
(3D) (solid dark-blue line) and Kolb (solid light-blue line) overpredict the 
measured v 2 - As these results do not include dissipative effects, such as those 
resulting from hadronic interactions in the final state, their failure indicates 
that further work will be necessary before a quantitative determination of the 
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viscosity in the QGP phase is possible. Progress will require both theoretical 
development and experimental measures that are less sensitive to how the 
azimuthal asymmetry of the energy-momentum tensor is distributed between 
different particles in the final stage of the reaction, e.g. the elliptic flow of the 
total transverse energy. 

The same hydro models are now compared to the measured spectra from cen- 
tral collisions. The bottom right panel of Fig. 20 shows that all the hydro 
models reproduce the pion spectra below p T ~ 1 GeV/c; at higher p T the 
particles are less likely to be equilibrated and hydro models are not expected 
to work well. In the bottom left panel the calculated proton spectra from 
Huovinen [107] (solid green line) are lower than the data, due to the calcu- 
lation maintaining chemical equilibrium throughout the hadronic phase. The 
lower temperature chemical freezeout suppresses the final calculated yield of 
heavier particles such as protons. Of the two calculations from Teaney [106] 
the calculation that includes the QGP phase (solid red line) reproduces the 
proton spectra, presumably because of the increased transverse flow from the 
stronger early pressure gradients. Hirano's and Kolb's (solid dark and light- 
blue lines) calculations break chemical equilibrium during the hadronic phase 
and overpredict the proton spectra at low px- 

One difficulty is that the spectra comparison with hydrodynamic models is 
for central collisions while the v 2 comparison is for minimum-bias collisions. 
It is dificult to use central collisions for the v 2 comparison since the collisions 
are nearly symmetric and hence v 2 is small. In addition, hydrodyamic calcu- 
lations that reproduce v 2 values over a broad range of centrality (from 0-45% 
in Ref. [103]) tend to overpredict the data for more peripheral collisions by 
approximately 25%, presumably because of a breakdown in the hydrodynamic 
assumptions. Hence when comparing to minimum-bias data sets, an overpre- 
diction of v 2 from the hydro models of less than 20% should be acceptable. 

These comparisons between data and hydro models are summarized in Table 
1 and in the following conclusions; 

• v 2 (pt, PID) is sensitive to all stages of the reaction. Elliptic flow is produced 
by strong scattering in the initial phase, while the detailed shape of v 2 (pr) 
and how the momentum asymmetry is distributed to different particles is 
affected by the transition from a QGP to hadronic phase and scattering in 
the final hadronic stage. 

• The hydro models that reproduce the low-pr proton v 2 are those that include 
both a QGP and hadronic phase. 

• The hadronic phase critically affects the final values of v 2 (px, PID). Models 
(Hirano, Kolb) that include partial chemical equilibrium to reproduce the 
baryon yield, completely fail on the pion v 2 . 

• The only model that survives this comparison with measured v 2 and spectra 



44 





QGP+mixed+RG 


mixed+RG 


RG 


Teaney 


Hirano 


Kolb 


Huovinen 


Teaney 


Huovinen 


Reference 


[1UDJ 




[1 1 7 1 931 


n D7i 

L- LU ' J 


[106] 


[107] 


latent 
heat 

(GeV/fm 3 ) 


0.8 


1.7 


1.15 


1.15 


U.o 




mit. Ernax 

(GeV/fm 3 ) 


16.7 




23 


23 


16.7 


23 


init. (e) 
(GeV/fm 3 ) 


11. U 


1 9 r 

13.5 






11.0 




To fm/c 


1.0 


0.6 


0.6 


U.D 


1.0 


0.6 


hadronic 
stage 


KL)1V1JJ 


partial 

chemical 

equil. 


partial 

chemical 

equil. 


lull equil. 


KLJMJJ 


lull equil. 


proton v2 


yes 


< 0.7 
GeV/c 


< 0.7 

GeV/c 


yes 


no 


no 


pion v2 


yes 


no 


no 


yes 


yes 


yes 


proton 
spectra 






l-J X K .L UU.1L Li 


L±\J 


no 


no 


pion spec- 
tra 


yes 


< 1 GeV/c 


< 1 GeV/c 


yes 


< 0.7 
GeV/c 


yes 


HBT 


Not avail- 
able 


No 


No 


No 


Not avail- 
able 


Not avail- 
able 



Table 1 



Summary of various hydro model assumptions and a comparison between measure- 
ments and hydro calculations. Two initial energies are tabulated, either the max- 
imum energy density at the center of the collision or the energy density averaged 
over the transverse profile. 

is Teaney's (solid red line) which includes a strong expansion in a QGP 
phase, a phase transition to a mixed phase, and then a hadronic cascade 
in the final hadronic state. There are open questions in this hybrid model, 
e.g. the sensitivity of the results to the matching conditions between hydro 
and RQMD. All other models fail in at least one v 2 or spectra comparison, 
partially due to differences in modeling the final hadronic state. 
• Until the model uncertainty in the final state is reduced, it is not yet possible 
to use the measured splitting between proton and pion ^(pt) to extract 
quantitative information on the EOS during the reaction, including the 
possible softening of the EOS due to the presence of a mixed phase. 
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A comparison with the HBT data and some of the hydro models is shown 
in Fig. 21. It is unfortunate that not all hydro models have been compared 
to HBT data, e.g. the hydro+RQMD model from Teaney [106] has not been 
confronted with this observable. The hydro calculation from Kolb, Heinz and 
Huovinen [123] (solid green line) includes a first-order phase transition which 
leads to a long lifetime for the system. The source parameter -Ri ong is considered 
most sensitive to the duration of the whole collision, i.e. from initial overlap to 
final particle emission, and the Kolb/Huovinen hydro calculation [123] (solid 
green line) overpredicts the measured -Ri ong data. Changing to partial chemical 
equilibrium in the hadronic stage [118], indicated with the dark blue line, 
reduces the lifetime of the collision which improves the agreement with -Ri on g- 
However the ratio Rout/ Rsidc, which is sensitive to the duration over which 
particles are emitted, is still overpredicted. 

There have been many attempts to understand what may be causing the 
disagreement with data (known collectively as the HBT puzzle) ; 

• Sinyukov et al. [124] and Grassi et al. [125] have suggested that the sharp 
Cooper-Frye freezeout condition [126] should be replaced by an emission 
function that decouples hadrons depending on their hadronic cross section. 

• However when this has been effectively implemented by using a hadronic 
cascade (URQMD) for the final hadronic stage, the predicted ratio R ou t/ R s \a c 
increases and diverges further from the data [127]. Modeling the final stage 
with a hadronic cascade effectively includes dissipative effects which should 
increase the duration of emission and produce a larger ratio -R out /-Rside • 

• One method to reduce the lifetime of the reaction is to change the QGP 
EOS. Using a crossover instead of a first-order transition reduces the ratio 
-Rout/-R S idc by about 50% to i? ou t/-Rside ~ 1-5 [128] which is still larger than 
the data. Because the calculation was restricted to rj = Zschiesche et al. 
were unable to compare with the measured values of -Ri ong - 

In summary, model comparisons seem to be closer to the HBT data when the 
lifetime of the collision is made smaller than the long time resulting from a 
first-order phase transition. The small values of i? ou t/-Rside may indicate that 
there is little to no mixed phase present in the reactions. One possible direction 
for future comparisons with data is to include a more realistic EOS into the 
hydro models, e.g. to take the EOS from lattice QCD calculations [129]. Such a 
calculation needs to be compared with all the available data, including spectra 
and v 2 , as well as HBT. 



3. 6 Conclusions 

In summary we can make the following conclusions 
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Fig. 21. The kx dependence of the Bertsch-Pratt parameters for tt + tt + (blue square) 
and tt~tt~ (red circle) for 0-30% centrality with statistical error bars and systematic 
error bands. Results from PHENIX [56], STAR [114] and hydrodynamics models 
(Hirano [118], Kolb/Huovinen [123] and Soff [130], diamonds) are overlaid. 



• The measured yields and spectra of hadrons are consistent with thermal 
emission from a strongly expanding source. 

• Strangeness is fully saturated at RHIC, consistent with full chemical equi- 
librium. 

• The scaling of v-i with eccentricity shows that collective behavior is estab- 
lished early in the collision. 

• Elliptic flow is stronger at RHIC than at the SPS, since the measured slope 
of v 2 (pt) f° r pions is 50% larger at RHIC. 

• The measured proton t^Pr) is less than that for pions at low pt] the small 
magnitude of the proton t> 2 at low px is reproduced by hydro models that 
include both a QGP and hadronic phase. 

• However several of the hydro models that reproduce the proton v 2 {pt) fail 
for the pion v 2 (pt)- 

• The HBT source parameters, especially the small value of -Ri ong and the ratio 
-R ou t/-R S ide, suggest that the mixed phase is too long-lived in the current 
hydro calculations. 
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Hence we currently do not have a consistent picture of the space-time dynamics 
of reactions at RHIC as revealed by spectra, t>2, and HBT. The lack of a 
consistent picture of the dynamics means that it is not yet possible to extract 
quantitative properties of the QGP or mixed phase using the observables v 2 , 
spectra, or HBT. 



4 FLUCTUATIONS 

4-1 Net- Charge Fluctuations 

In the study of the fluctuations of multiplicity as a means to understand the 
dynamics of charged particle production, one important realization was to 
use small regions of phase space, where energy-momentum conservation con- 
straints would not be significant [131,132,133]. Such studies led to the impor- 
tant observation that the distribution of multiplicity, even in small intervals 
near mid-rapidity, was Negative Binomial rather than Poisson, which indi- 
cated large multiplicity correlations even in small 8r] intervals [134]. No such 
studies are yet available at RHIC. 

Based on predictions that event-by-event fluctuations of the net charge in 
local phase space regions would show a large decrease as a signature of the 
QGP [135,136,137], net-charge fluctuations were measured in PHENIX [44]. 
The idea is that in a QGP composed of fractionally charged quarks, the larger 
number of fractionally charged particles compared to unit-charged hadrons 
would result in smaller relative net-charge fluctuations in a QGP than for 
a pure gas of hadrons and that this original fluctuation would survive the 
transition back to ordinary hadrons. 

It is important to realize that the study of net-charge fluctuations represents 
the study of fluctuations in a quantity that is conserved over all phase space. 
Consider N = N + + iV_ charged particles produced in the full phase space. 
By charge conservation N + = AL = N/2, and the net charge Q = N + — AL 
is identically zero so that there are no net-charge fluctuations — the variance 
V(Q) = 0, where 

V(Q) = (Q 2 ) - (Q) 2 . (10) 

In a smaller region of phase space, where p is the fraction of iV observed in a 
stochastic scenario, the mean and variance of the number of positive n+ and 
negative n_ particles are equal, but the variance of Q is no longer identically 
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zero: 



{n+) = {n-)=pN/2, 



(11) 



V(n + ) = V(n-)=p(l-p)N/2, 



(12) 



from which it follows that 



V(Q) = (l-p)n ch , 



(13) 



where n ch = pN is the expected number of charged particles on the interval. 
Thus the normalized variance in Q (normalized to Poisson statistics) is defined 
as: 



In the limit n c h 0, the variance of the charge ratio R — n + /n- approaches 
V(R) = 4(l—p)/n ch . However, it is well known in mathematical statistics that 
moments of the inverse of a stochastic variable, e.g. l/n_, diverge if there is 
any finite probability, no matter how small, for n_ = 0. Thus, the charge ratio 
is not a stable measure of fluctuations. 

The previous arguments are based on fixed N. The results where N varies 
according to a specified distribution are also interesting. If n_ is Poisson dis- 
tributed, with mean value \x = N/2 over the whole phase space, then in the 
region of phase space with probability p the distribution is also Poisson, with 
mean = ji p = pN/2. If, on the other hand n_ is Negative Binomial dis- 

tributed, with mean value ix — N/2 and NBD parameter a 2 / fi 2 — l/ii = 1/k 
for the whole phase space, then in the region of phase space with probability 
p, the distribution is Negative Binomial with mean (n-)\ p — \x v — pN/2 and 
the same value of 1/k. 

Actually, the binomial division preserves <j 2 /fi 2 — l/fi p = 1/k, for any distribu- 
tion [138]. This appears to indicate that smaller intervals, which tend to have 
larger values of <J 2 /n p would be less sensitive to the global 1/k, the long-range 
correlation. This would be true except for the fact that there are short-range 
correlations which are better seen on small intervals of phase space. Another 
important thing to note regarding a binomial split of a Negative Binomial 
distribution is that the two subintervals are not statistically independent. The 
conditional probability distribution on the interval (1 — p) depends upon the 
outcome on the interval p [139]. It is unfortunate that these elegant arguments 
can not be applied to the net-charge fluctuations since (Q) = 0. 



v(Q) 



V(Q) 



(l-p). 



(14) 



n c h 
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The PHENIX measurement [44] of the normalized variance v(Q) of net-charge 
fluctuations is shown in Fig. 22 in the interval —0.35 < 5i] < +0.35 as a 
function of the azimuthal angular interval of reconstructed tracks, either at 
the detector, A(p d , or at the vertex, A0 r , chosen symmetrically around the 
detector acceptance. For smaller A0 r the data agree with the purely stochastic 
(1 — p) dependence shown as the solid line, but deviate from the stochastic 
prediction at larger values due to correlations from resonance decay, such as 
p° -> 7T+ + 7r~ as nicely explained by RQMD [140]. 

Absent new theoretical insight, it is difficult to understand how quark- level 
net-charge fluctuations in a QGP can be related to net-charge fluctuations of 
hadrons, where, by definition, strong correlations exist, e.g., in the formation 
of a meson from a q — q pair. Also, the study of the fluctuations of net charge, 
which is conserved, may not be as useful to detect interesting fluctuations as 
the study of fluctuations of the total charged multiplicity, which is much less 
constrained by conservation laws. This has yet to be tried at RHIC. 

4-2 Event-by-Event Average-pr Fluctuations 

Fluctuations in the event-by-event average pr, denoted M PT , have been mea- 
sured and provide a severely small limit on possible fluctuations from a sharp 
phase transition. For events with n detected charged particles with magnitudes 
of transverse momenta, px v the event-by-event average px, denoted M PT is de- 
fined as: 

1 n 

M PT =pf=-J2PT t . (15) 
n i=i 

Mixed events are used to define the baseline for random fluctuations of M PT in 
PHENIX [45,55]. This has the advantage of effectively removing any residual 
detector-dependent effects. The event-by-event average distributions are very 
sensitive to the number of tracks in the event (denoted n or N trac k s ), so the 
mixed event sample is produced with the identical N tra cks distribution as the 
data. Additionally, no two tracks from the same data event are placed in 
the same mixed event in order to remove any intra-event correlations in px- 
Finally, (M PT ) must exactly match the semi-inclusive (pt)- 

For the case of statistical independent emission, where the fluctuations are 
purely random, an analytical formula for the distribution in M pT can be ob- 
tained assuming Negative Binomial distributed event-by-event multiplicity, 
with Gamma distributed semi- inclusive pt spectra [141]. The formula depends 
on the four semi-inclusive parameters (n), 1/k, b and p which are derived from 
the means and standard deviations of the semi-inclusive pt and multiplicity 
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Fig. 22. v(Q) for the 10% most central events in data and RQMD, as a function of 
the azimuthal interval in degrees of reconstructed tracks, either a) at the detector, 
A(f>d, or b) at the vertex, A(j) r , chosen symmetrically around the detector acceptance. 
For data, the error band shows the total statistical error, whereas the error bars 
indicate the uncorrelated part. The solid line shows the expected reduction in v(Q) 
in the stochastic scenario when global charge conservation is taken into account 
[44]. 
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distributions, (n), a n , (px), & PT ' 

^max 

f(y)= fNBi>(n,l/k,(n))fr(y,np,nh), (16) 



where y = M PT . For fixed n, and purely random fluctuations, the mean 
and standard deviation of M PT follow the expected behavior, (M PT ) = (px), 
a M PT = °p T /v^- I 11 PHENIX, Eq. 16 is used to confirm the randomness of 
mixed events which are used to define the baseline for random fluctuations of 
M PT [45,55]. 

The measured M PT distributions for the data in two centrality classes for 
v /s^7 = 200 GeV Au+Au collisions in PHENIX [55] are shown in Fig. 23 
(data points) compared to the mixed-event distributions (histograms). The 
non-Gaussian, Gamma-distribution shape of the M PT distributions is evident. 
The difference between the data and the mixed-event random baseline distri- 
butions is not visible to the naked eye. The nonrandom fluctuation is quantified 
by the fractional difference of u PT , the normalized standard deviation of M PT , 
for the data and the mixed-event (random) samples: 

t~t ^Pt, data ^Pt, mixed /io\ 

t PT = - . (18) 

^Pt, mixed 



The results are shown as a function of centrality, represented by N part in Fig. 
24. 

The dependence of F PT on N part is striking. To further understand this de- 
pendence and the source of these nonrandom fluctuations, F PT was measured 
over a varying px range, 0.2 GeV/c < px < p™ ax (Fig- 25), where p™ ax = 2.0 
GeV/c for the N part dependence. 

The increase of F PT with p™ ax suggests elliptic flow or jet origin. This was 
investigated using a Monte Carlo simulation of correlations due to elliptic 
flow and jets in the PHENIX acceptance. The flow was significant only in 
the lowest centrality bin and negligible (F PT < 0.1%) at higher centralities. 
Jets were simulated by embedding (at a uniform rate per generated particle, 
S P rob(N part )) p + p hard-scattering events from the PYTHIA event generator 
into simulated Au+Au events assembled at random according to the measured 
Ntracks and semi-inclusive px distributions. This changed (px) and a PT by less 
than 0.1%. S pro b(N part ) was either constant for all centrality classes, or scaled 
by the measured hard-scattering suppression factor RAA(N part ) for px > 4.5 
GeV/c [49]. A value F PT = 2.06% for p + p collisions was extracted from 
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Fig. 23. M pr for 30-35% and 0-5% centrality classes: data (points) mixed events 
(histogram) [55]. 
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Fig. 24. F PT vs. centrality, represented as the average number of participants (N part ) 
in a centrality class, compared to jet simulation [55]. 

pure PYTHIA events in the PHENIX acceptance in agreement with the p + 
p measurement (Fig. 24). The value of S pro b(N part ) was chosen so that the 
simulation with S pro b(N part ) x RAA(N part ) agreed with the data at N part = 
182. The centrality and p™ ax dependences of the measured F PT match the 
simulation very well, but only when the Raa scaling is included. 

A less experiment-dependent method to compare nonrandom fluctuations is 
to assume that the entire F PT is due to temperature fluctuations of the initial 
state, with RMS variation ctt/(T) [142,45]. Then, 



2 — 2 — fi — 

^PTjdata ^Pti mixed \^ 



a. 



r 



n)'{TY 



^Fpr^pT .mixed > 



(19) 



This yields cry/(T) = 1.8% for central collisions in PHENIX with similarly 
small values for the other Relativistic Heavy Ion experiments [143], 1.7% in 
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Fig. 25. F Pt vs. p™ ax compared to iZ^A-scaled jet simulation for the 20-25% cen- 
trality class (N part = 181.6) [55]. 

STAR, 1.3% in CERES, and 0.6% in NA49. These results put severely small 
limits on the critical fluctuations that were expected for a sharp phase tran- 
sition, both at SPS energies and at RHIC, but are consistent with the expec- 
tation from lattice QCD that the transition is a smooth crossover [12]. 

Other proposed explanations of the centrality and p™ ax dependences of F PT 
include: overlapping color strings which form clusters so that the number of 
sources and (pt) per source is modified as a function of centrality [144]; and 
near equilibrium p T correlations induced by spatial inhomogeniety [145]. 



4-3 Conclusions 



Critical behavior near the phase boundary can produce nonrandom fluctua- 
tions in observables such as the net-charge distribution and the average trans- 
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verse momentum distribution. Our search for net-charge fluctuations has ruled 
out the most naive model of charge fluctuations in a QGP, but it is unclear 
whether the charge fluctuation signature can survive hadronization. Our mea- 
surement of the event-by-event average pt distribution shows a nonrandom 
fluctuation that is consistent with the effect expected from high-p T jets. This 
puts a severe constraint on the critical fluctuations that were expected for 
a sharp phase transition but is consistent with the expectation from lattice 
QCD that the transition is a smooth crossover [12]. 



5 BINARY SCALING 



5.1 Hard Scattering and pQCD 



One way to get a partonic probe into the midst of an A + A collision is to use 
the high-pT partons produced by hard scattering. For p + p collisions in the 
RHIC energy range, hard scattering is considered to be the dominant process 
of particle production with px > 2 GeV/c at mid-rapidity. Typically, particles 
with px > 2 GeV/c are produced from states with two roughly back-to-back 
jets which are the result of scattering of constituents of the nucleons (partons) 
as described by pQCD [146]. 

The overall p+p hard-scattering cross section in "leading logarithm" pQCD is 
the sum over parton reactions a+b — > c+d (e.g. g+q — > g+q) at parton-parton 
center-of-mass (cm.) energy V^: 



d 3 a _ ly, , ttc^(Q 2 
dxidx 2 dcos9* s ^ 1x\X 2 



E/a(^l)/6(^)^^S ab (c0Sr) , (20) 



where f a (xi), fb(x 2 ), are parton distribution functions, the differential proba- 
bilities for partons a and b to carry momentum fractions X\ and x 2 of their re- 
spective protons (e.g. u(x 2 )), and where 9* is the scattering angle in the parton- 
parton cm. system. The parton-parton cm. energy squared is s = X\x 2 s, 
where y/s is the cm. energy of the p + p collision. The parton-parton cm. 
system moves with rapidity y = l/2\n(x 1 /x 2 ) in the p + p cm. system. 

Equation 20 gives the pr spectrum of outgoing parton c, which then fragments 
into hadrons, e.g. 7T°. The fragmentation function (z, fi 2 ) is the probability 
for a 7T° to carry a fraction z = p w /p c of the momentum of outgoing parton c. 
Equation 20 must be summed over all subprocesses leading to a 7r° in the final 
state. The parameter fi 2 is an unphysical "factorization" scale introduced to 
account for collinear singularities in the structure and fragmentation functions 
[146,147]. 
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In this formulation, / a (xi,/i 2 ), fb(x2,fi 2 ) and D^(z,fi 2 ) represent the "long- 
distance phenomena" to be determined by experiment; while the character- 
istic subprocess angular distributions, X afc (cos#*), and the coupling constant, 
a s (Q 2 ) = ^ln(Q 2 /A 2 ), are fundamental predictions of QCD [148,149,150] 
for the short-distance, large-Q 2 , phenomena. The momentum scale Q 2 ~ p\ 
for the scattering subprocess, while Q 2 ~ s for a Compton or annihilation 
subprocess, but the exact meaning of Q 2 tends to be treated as a parameter 
rather than a dynamical quantity. The transverse momentum of a scattered 
constituent is: 

Vt = Vt = \ sine*. (21) 



Equation 20 leads to a general '^-scaling' form for the invariant cross section 
of high-p T particle production: 



where x T = 2p T /y/s. The cross section has two factors, a function F(x T ) 
(G(xt)) which 'scales', i.e. depends only on the ratio of momenta, and a di- 
mensioned factor, 1/pt (l/y^"), where n equals 4 in lowest-order (LO) cal- 
culations, analogous to the l/q A form of Rutherford Scattering in QED. The 
structure and fragmentation functions are all in the F(xt) (G(xt)) term. Due 
to higher-order effects such as the running of the coupling constant, a s (Q 2 ), 
the evolution of the structure and fragmentation functions, and the initial- 
state transverse momentum kr, n is not a constant but is a function of xt, 
y/s. Measured values of ti(xt, y/s) in p + p collisions are between 5 and 8. 



5.2 Mid- Rapidity p T Spectra from p + p Collisions 



The scaling and power-law behavior of hard scattering are evident from the 
y/s dependence of the pt dependence of the p+p invariant cross sections. This 
is shown for nonidentified charged hadrons, (h + + h~)/2, in Fig. 26a. At low 
Pt < 1 GeV/c the cross sections exhibit a "thermal" exp (—Qpr) dependence, 
which is largely independent of y/s, while at high pt there is a power-law 
tail, due to hard scattering, which depends strongly on y/s. The characteristic 
variation with y/s at high p T is produced by the fundamental power-law and 
scaling dependence of Eqs. 20, 22. This is best illustrated by a plot of 

^n^v-s) x E ^a = QM ^ (23) 
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Fig. 26. (left) Ed?a(pT)/d 3 p at mid-rapidity as a function of yfs in p + p collisions, 
(right) y^GeV) 6 ' 3 x Ed?a/d?p vs. xt = 2px/\/s [53] (and references therein). 

as a function of xt, with n(xT, y/s) = 6.3, which is valid for the xt range 
of the present RHIC measurements (Fig. 26b). The data show an asymptotic 
power law with increasing xt- Data at a given yfs fall below the asymptote 
at successively lower values of xt with increasing y/s, corresponding to the 
transition region from hard to soft physics in the px region of about 2 GeV/c. 

The PHENIX measurement of the invariant cross section for 7T° production 
in p + p collisions at y/s = 200 GeV [60] agrees with NLO pQCD predictions 
over the range 2.0 < pr < 15 GeV/c (Fig. 27). 



5. 3 Scaling Hard Scattering from p + p to p + A and A + B Collisions 



Since hard scattering is point like, with distance scale 1/pr < 0.1 fm, and the 
hard-scattering cross section factorizes as shown in Eq. 20, the cross section 
in p + A or A + B collisions, compared to p + p, is proportional to the relative 
number of possible point-like encounters. The number of encounters of point- 
like constituents of nucleons is then proportional to A (AB), for p + A (A + B) 
minimum-bias collisions. For A + B collisions at impact parameter b, it is 
proportional to TabQ>), the nuclear thickness function, which is the integral 
of the product of nuclear thickness over the geometrical overlap region of the 
two nuclei. In detail, the semi-inclusive invariant yield of e.g. high-p^ 7r°'s for 
A + B inelastic collisions, with centrality /, is related to the p+p cross section 
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Fig. 27. PHENIX vr° invariant cross section at mid-rapidity from p + p collisions at 
y/s = 200 GeV, together with NLO pQCD predictions from Vogelsang [151,152]. 
a) The invariant differential cross section for inclusive tt° production (points) and 
the results from NLO pQCD calculations with equal renormalization and factor- 
ization scales of pt using the "Kniehl-Kramer-Potter" (solid line) and "Kretzer" 
(dashed line) sets of fragmentation functions, b) The relative statistical (points) 
and point-to-point systematic (band) errors. c,d) The relative difference between 
the data and the theory using KKP (c) and Kretzer (d) fragmentation functions 
with scales of j>t/2 (lower curve), pr, and 2px (upper curve). In all figures, the 
normalization error of 9.6% is not showgc[60]. 



by: 



d 2 N£ B 



N% dp T dy 



Note that 



Irp v JjTABiVfb _ ( Ncoll ) f 

{lAB)f ~ f f (l-e-°™ T AB{b))d 2 b - aNN > W 



where (N co u) f is the average number of binary nucleon-nucleon inelastic col- 
lisions, with cross section oattv, m the centrality class /. This leads to the 
description of the scaling for point-like processes as binary-collision (or N co u) 
scaling. 

Nuclear medium effects, either in the initial or final state, can modify the 
expected scaling. These modifications can be quantitatively studied by mea- 
surement of the Nuclear modification factor Rab, which is defined as 

R = dN ^ = dN %* (2(]] 

{Tab) f x da p NN (N coll ) f x dNg N ' 1 ° J 



where dN^ B is the differential yield of a point-like process P in a A+B collision 
and da^N * s the cross section of P in N + N collision. If there are no initial- 
or final-state effects that modify the yield of P in A + B collisions, the process 
P scales with (Tab) / an d Rab = 1- Sometimes, the central to peripheral 
ratio, Rcp, is used as an alternative to Rab- The central to peripheral ratio 
is defined as 



^]\j~Central j ^jy Central^ 
^Peripheral j ^Peripheral ^ 



/ \ coil / (2,7) 

CP i atd™^!,™/ 1 1 njPeripheral\ ' v / 



where dN Central and dJSf Pen vherai are ^ e differential yield per event of the 
studied process in a central and peripheral collision, respectively. If the yield 
of the process scales with the number of binary collisions, Rcp = 1. 



5.4 Binary Scaling in I + A, p + A, and Low-Energy A + A 



In deeply inelastic lepton scattering, where hard scattering was discovered 
[153,154,155], the cross section for /i-A collisions is indeed proportional to 
A im (Fig. 28). This indicates that the structure function of a nucleus of mass 
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Fig. 28. [i- A cross section vs. A [156]. 

A is simply A times the structure function of a nucleon (with only minor 
deviations, < 10% for 0.02 < x < 0.50 [157]), which means that the nucleus 
acts like an incoherent superposition of nucleons for hard scattering of leptons. 

The situation is rather different in p + A collisions: the cross section at a 
given p T also scales as a power law, A 01 ^^ (Fig. 29), but the power a(px) 
is greater than 1. This is called the "Cronin Effect" [158]. The enhancement 
(relative to A im ) is thought to be due to the multiple scattering of the incident 
partons while passing through the nucleus A before the collision [159,160], 
which smears the axis of the hard scattering relative to the axis of the incident 
beam, leading to the characteristic "Cronin Effect" shape for R A (p T ) (Fig. 
30). At low pt < 1 GeV/c, the cross-section is no longer point like, so the 
scattering is shadowed (oc A 2 ^ 3 ), thus Ra < 1. At larger p T > 2 GeV/c, as 
the hard-scattering, power-law p T spectrum begins to dominate, the multiple 
scattering smears the spectrum to larger p T leading to an enhancement relative 
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Fig. 29. Cronin Effect mp + A, for tt~ with p T = 4.61 GeV/c. a(pr) = 1.148±0.010 
[158]. 

to binary-scaling which dissipates with increasing px as the influence of the 
multiple scattering diminishes. 

Previous measurements of high-p^ particle production in A + A collisions at 
^/s NN < 31 GeV (Fig. 31) and in p + A (or d + A) collisions (Fig. 30) including 
measurements at RHIC [64] at mid-rapidity (Fig. 32) all show binary scaling 
or a Cronin effect. 

This establishes that the initial condition for hard scattering at RHIC at 
mid-rapidity is an incoherent superposition of nucleon structure functions, 
including gluons, where multiple scattering before the hard collision smears 
the pt spectrum of scattered particles to be somewhat above the simple point- 
like binary (N co u) scaling. 
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Fig. 30. Cronin effect at fixed target energies expressed as -R\y/Be> the ratio of the 
point-like scaled cross sections in p + W and p + Be collisions vs. pr [161]. 

An alternative view of the initial state of a nucleus at RHIC is provided by 
the color glass condensate (CGC), in which the gluon population at low x is 
not an incoherent superposition of nucleon structure functions but is limited 
with increasing A by non-linear gluon-gluon fusion resulting from the overlap 
of gluons from several nucleons in the plane of the nucleus transverse to the 
collision axis [168]. A Cronin effect in d+A collisions, as shown in Fig. 32, can 
be reproduced in the CGC with a suitable choice of initial state parameters, 
which must also reproduce quantitatively the observed binary scaling of the 
direct photon production and total charm production in Au+Au collisions to 
be shown below (Figs. 33, 34). However, at this writing, no detailed quantita- 
tive description of the CGC initial state which satisfies these three conditions 
has been published. 



5.5 Binary Scaling in Au+Au Collisions at RHIC — Direct Photons and Charm 
Yield 



The production of hard photons in Au+Au collisions at RHIC via the con- 
stituent reactions (e.g. g + q — > 7 + q) is a very important test of QCD and 
the initial state, because the photons only interact electromagnetically, hence 
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Fig. 31. Nuclear modification factors for ir° production at the CERN-ISR in min- 
imum-bias a + a reactions at ■ v /ijv7V = 31 GeV [162] and for pion production at 
the CERN-SPS in central Pb+Pb [163], Pb+Au [164], and S+Au [165] reactions at 
\J s nn ~ 20 GeV. The Raa from SPS are obtained using the p + p parametrization 
proposed in ref. [166]. The shaded band around Raa = 1 represents the overall frac- 
tional uncertainty of the SPS data (including in quadrature the 25% uncertainty of 
the p + p reference and the 10% error of the Glauber calculation of N co u). There 
is an additional overall uncertainty of ±15% for the CERES data not shown in the 
plot [164]. 



hardly at all, with any final-state medium produced. The direct-photon cross 
section and centrality dependence should then reflect only the properties of the 
initial state, notably the product of the gluon and quark structure functions 
of the Au nuclei. 



The first measurement of direct photon production in Au+Au collisions at 
RHIC has been reported by the PHENIX collaboration (Fig. 33) [73]. The 
data exhibit pure point-like (N co u) scaling as a function of centrality relative 
to a pQCD calculation for p + p collisions. The statistical and systematic 
errors still leave some room for a small Cronin effect and/or some thermal 
photon production. The observation of direct photon production establishes 
the importance of gluon degrees of freedom at RHIC. 

PHENIX measured the single-electron yield from nonphotonic sources in Au+Au 
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Fig. 32. Cronin effect in Rcp, the ratio of point-like scaled central to peripheral 
collisions for pions in d+Au at y / s/v r /v = 200 GeV[167]. Data points for low px are 
7r identified by Time of Flight (TOF). Data at medium px are for ir° identified by 
reconstruction in the Electromagnetic Calorimeter (EMCAL). Highest pt data are 
for ir^ identified by a count in the Ring Imaging Cerenkov Counter (RICH) and a 
deposited energy /momentum and shower shape in the EMCAL inconsistent with 
those of a photon or electron. The shaded band on the right represents the overall 
fractional systematic uncertainty due to N co u. 



collision at 130 GeV [42] and 200 GeV [57]. Since semi-leptonic decay of 
charm is the dominant source of the non-photonic electrons at low px (pt < 3 
GeV/c), the total yield of charm can be determined from the integrated yield 
of non-photonic electrons in the low-px region. Figure 34 shows the yield of 
non-photonic electrons (0.8 < px < 4.0 GeV/c) per NN collision in Au+Au 
reactions at ^snn = 200 GeV as a function of N co u [57]. The N co u dependence 
of the yield is fit to N^ oll , where a = 1 is the expectation for binary scaling. 
We find a = 0.938 ± 0.075(stat.)±0.018(sys.), showing that the total yield of 
charm-decay electrons is consistent with binary scaling. It should be noted 
that medium effects, such as energy loss of charm in the dense hot medium, 
can only influence the momentum distribution of charm, and have little effect 
on the total yield of charm. Initial-state effects, such as shadowing, and other 
effects, such as thermal production of charm, are believed to be very small 
for charm production at RHIC energy. Therefore, the observation of binary 
scaling of the total charm yield in Au+Au collisions at RHIC may also be 
considered as an experimental verification of the binary scaling of a point-like 
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Fig. 33. PHENIX direct photon measurements relative to the background for for 
minimum bias and for five centralities of Au+Au collisions at ^/s NN = 200 GeV 
(0-10% is the most central). Statistical and total errors are indicated separately 
on each data point by the vertical bar and shaded region, respectively. The curves 
represent a pQCD calculation of direct photons in p + p collisions from Vogelsang 
[169,170,171,172] scaled to Au+Au assuming pure point-like (iV co ii) scaling, with 
no suppression. The shaded region around the curves indicate the variation of the 
pQCD calculation for scale changes from pr/2 to 2px, plus the (N co u) uncertainty 
[73]. 

pQCD process. 



5. 6 Conclusions 



In this section evidence has been presented to show that the initial condition 
for hard-scattering at RHIC at mid-rapidity is an incoherent superposition of 
nucleon structure functions, including gluons, where multiple scattering be- 
fore the hard collision can smear the px spectrum of scattered particles to be 
somewhat above the simple point-like binary (N co u) scaling. This was demon- 
strated using the reactions: pion production in d+Au collisions, where there 
is no final-state medium, and direct photon production in Au+Au collisions, 
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Fig. 34. Non-photonic electron yield (0.8 < px < 4.0 GeV/c), dominated by 
semi-leptonic charm decays, measured in Au+Au collisions at ^/sJm = 200 GeV 
scaled by N co u as a function of N co u. The right-hand scale shows the corresponding 
electron cross section per NN collision in the above pt range. The yield in p + p 
collision at 200 GeV is also shown [57]. 

where the outgoing photons interact electromagnetically, hence hardly at all, 
with any final-state medium produced. The total charm yield in Au+Au, a 
reaction dominated by the subprocess g + g — > c + c, and which is not sensitive 
to final-state medium effects for the total yield of c + c pairs, also exhibits 
binary scaling. The latter two measurements provide experimental evidence 
for the binary scaling of point-like pQCD processes in Au+Au collisions. 

The color glass condensate (CGC) provides an alternative view of the initial 
state of a nucleus at RHIC in which coherence of gluons due to non-linear 
gluon-gluon fusion can produce a Cronin-like effect, depending on the initial 
conditions and the kinematic range covered. However, at the present writing, 
there is no CGC description of the initial state nuclear structure function 
which reproduces the observed Cronin effect for pions in d+Au collisions and 
the observed binary scaling for both direct photon production and the total 
charm yield in Au+Au collisions. 



6 NUMBER DENSITY AND HIGH p T SUPPRESSION 



To study the initial properties of the matter created in heavy ion collisions 
we need a probe that is already present at earliest times and that is directly 
sensitive to the properties of the medium. Partons resulting from hard scat- 
terings during the initial crossing of the two nuclei in A+A collisions provide 
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such a probe. Energetic partons propagating through a dense medium are 
predicted to lose energy [173,174,175,176,177,178,179,180,181] thus producing 
a suppression in the yield of high-p^ hadrons produced from the fragmen- 
tation of these partons. Initial measurements from RHIC Run 1 [39,47,182] 
and Run 2 [49,53,183,184] demonstrated such a suppression, and the results 
of d+Au measurements [64,68,67,66] showed that the suppression was not due 
to initial-state effects. Further measurements have indicated a modification 
of di-jet angular correlations [68] that has also been attributed to in-medium 
parton energy loss [185,186]. 

While the energy loss of hard-scattered partons was originally proposed as a 
signature of the quark-gluon plasma and deconfinement, it has been argued 
recently that the energy loss is sensitive only to the density of unscreened color 
charges and not directly to deconfinement [177,178,179,180,181,187,188]. Ide- 
ally, a measurement of initial parton densities together with constraints on ini- 
tial energy density might allow an estimate of the temperature of the medium. 
As will be seen below, the current high-p^ measurements and theoretical tools 
for interpreting the experimental data are not yet sufficient to take such a step. 
Instead, the energy loss results are currently being used to provide estimates of 
the initial energy density. The remainder of this section summarizes PHENIX 
experimental data related to high-p T suppression, discusses the current state 
of theoretical understanding of the energy loss process and concludes with a 
statement of estimates for initial parton number and energy densities that 
currently can be made. 



6. 1 Single particle spectra, Raa 

As described in Section 5, in the absence of modifications due to initial-state 
or final-state effects, the rate for the production of particles through hard- 
scattering processes in nucleus-nucleus collisions is expected to be given by 
the equivalent p + p hard-scattering cross section multiplied by Tab- Figure 
35 shows PHENIX n° spectra, d 2 N/dpTdy, measured in 200 GeV [49] periph- 
eral (80-92%) and central (0-10%) Au+Au collisions compared to measured 
[60] p + p cross sections multiplied by the peripheral and central Tab values 
estimated using the procedure described in Section 5. The error bands on the 
p + p data points reflect both the systematic errors on the p + p cross sections 
and the uncertainties in the Tab values. As the figure clearly demonstrates, the 
central Au+Au n° yields are strongly suppressed relative to the "expected" 
yields over the entire measured pr range. In contrast, the peripheral yields 
compared to the T^s-scaled p + p cross sections show little or no suppression. 
The results incontrovertibly demonstrate that there is a strong and centrality- 
dependent suppression of the production of high-pr pions relative to pQCD- 
motivated expectations. This is quite different from measurements of Raa 
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in Pb+Pb collisions at y/s^N = 17.3 GeV where in semi-peripheral Pb+Pb 
collisions there is a nuclear enhancement increasing with px similar to the 
well-known Cronin effect, while in central collisions the Cronin enhancement 




p T (GeV/c) 

Fig. 35. 7T° pt spectra in 200 GeV Au+Au collisions [49] compared to a Tab scaling 
of the 200 GeV p + p 7r° differential cross section [60] . The central data were obtained 
with a 0-10% centrality cut while the peripheral data were obtained with an 80-92% 
cut. 

To better demonstrate quantitatively the suppression in central collisions in- 
dicated in Fig. 35, we show in Fig. 36 Raa(pt) for mid-rapidity 7r°'s in central 
and peripheral 200 GeV Au+Au collisions. We also show the values obtained 
from minimum-bias 200 GeV <i+Au collisions [64] which provide a stringent 
test of the possible contribution of initial-state nuclear effects to the observed 
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suppression in Au+Au collisions. The error bands on the data indicate com- 
bined statistical and point-to-point systematic errors and the bars shown next 
to the different data sets indicate common systematic errors due to uncertain- 
ties in the p + p cross section normalization and T AB . 

Figure 36 shows that the central Au+Au 7r° suppression changes only slightly 
over the measured pt range and reaches an approximately pr-independent 
factor of 5 (Raa ~ 0.2) for p T > 4 — 5 GeV/c. The peripheral Au+Au Raa 
values are consistent with one after taking into account systematic errors but 
we cannot rule out a slight suppression suggested by the peripheral Raa val- 
ues. In all of the data sets Raa decreases with decreasing pt for px < 2 GeV/ c. 
This decrease, known since the original measurements of the A dependence of 
particle production in pA collisions is due to contributions of soft hadronic pro- 
cesses at low p T that are expected to increase more slowly than proportional 
to Tab- The d+Au RdA values are also consistent with one within systematic 
uncertainties, but in contrast to the Au+Au results, the data suggest a slight 
enhancement. The ci+Au RdA values above 2 GeV/c exceed one for nearly the 
entire experimentally covered px range. As shown previously in Fig. 32, only 
for Pt^G GeV/c does the ci+Au pion yield return to the T^-scaling expecta- 
tion. Such a small enhancement is consistent with expectations based on prior 
measurements of the Cronin effect [189,190], and it is also quantitatively con- 
sistent with calculations incorporating the initial-state multiple scattering that 
is thought to produce the Cronin effect [191,192,193,194,195,196,197]. There- 
fore the Cronin effect at RHIC cannot mask a strong initial-state suppression 
of the parton distributions in the Au nucleus [198]. 

To better demonstrate the systematic behavior of the high-pT suppression 
we show in Fig. 37 7r° [49] and unidentified charged particle Raa values [53] 
as a function of Pt for various centrality bins. While for moderate px val- 
ues (2 < p T < 5 GeV/c) total charged particle production is suppressed less 
than pion production, the charged particle and n° Raa values become equal, 
within errors, at high px- This evolution in the charged particle suppression 
is related to contributions from the (anti)protons that will be discussed fur- 
ther below. Despite the differences resulting from the protons, the charged 
particles and 7r°'s exhibit very similar trends in the suppression vs. px and 
vs. centrality. The suppression increases smoothly with centrality though the 
change in Raa values at high px is most rapid in the middle of the centrality 
range. Figure 37 also shows that the suppression is approximately constant as 
a function of px for p T > 4.5 GeV/c in all centrality bins. We take advantage 
of this feature of the data to better illustrate the centrality dependence of the 
suppression by integrating both the Au+Au spectra and the reference p + p 
cross sections over pt > 4.5 GeV/c and using these integrated quantities to 
determine an average suppression factor, R A a for px > 4.5 GeV/c. We plot 
the charged particle and n° Raa values vs. N part in Fig. 38(top). This figure 
suggests that the suppression evolves smoothly with N pait , showing no abrupt 
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Fig. 36. 7T° Raa(pt) for central (0-10 %) and peripheral (80-92 %) Au+Au collisions 
[49] and minimum-bias d+Au collisions [64] . The shaded boxes on the left show the 
systematic errors for the Au+Au Raa values resulting from overall normalization 
of spectra and uncertainties in Tab- The shaded box on the right shows the same 
systematic error for the <i+Au points. 

onset of suppression. The charged particles and 7r°'s exhibit similar evolution 
of suppression with N part . In the most central collisions we obtain Raa val- 
ues of 0.24 ± 0M(total) and 0.23 ± 0.05(total) for charged particles and 7r 0, s 
respectively. In peripheral collisions, Raa approaches one, but the systematic 
errors on the most peripheral Tab values are sufficiently large that we cannot 
rule out ~ 20% deviations of the peripheral Au+Au hard-scattering yields 
from the T^s-scaled p + p cross sections. 



An alternative method for evaluating the evolution of the high-pr suppression 
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Fig. 37. Centrality and pr dependence of nuclear modification factors in 200 GeV 
Au+Au collisions [53]. Top panel: rr and charged particle Raa(pt) for ten centrality 
bins. Bottom panel: charged particle Rcp vs. pt- 
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Fig. 38. Top panel: Raa vs. iV par t obtained from p^-integrated (px > 4.5 GeV/c) 
Au+Au 7T° and charged-hadron spectra. The band indicates the systematic error 
bands on a hypothetical Tab scaling of the p + p pr-integrated cross section. Bottom 
panel: 7r° and charged hadron yield per participant vs. -/V par t divided by the same 
quantity in p + p collisions (R AA art ). The solid band shows the same band as in the 
top panel expressed in terms of yield per participant pair while the dashed band 
indicates the systematic error bands around a hypothetical iVp ar t scaling. Both plots 
are from [53]. 
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with centrality is provided in Fig. 38 (bottom) which presents the charged and 
7T° yields per participant integrated over p? > 4.5 GeV/c as a function of 
iV pa rt [53] divided by the same quantity in p + p collisions. Also shown in the 
figure are curves demonstrating the iVp art dependence that would result if the 
7T° and charged particle yields exactly Tab scaled and what an iV pa rt scaling 
from p + p collisions would imply. As Fig. 38 demonstrates, the high-pr yields 
of both charged hadrons and 7r°'s per participant increase proportional to Tab 
for small A^ part but level off and then decrease with increasing iV part in more 
central collisions. The PHENIX measurements do not naturally support an ap- 
proximate iV part scaling of high-p T particle production suggested in an analysis 
of PHOBOS data. The PHENIX R^a* values decrease from mid-peripheral 
(iVp ar t ~ 75) to central collisions by an amount larger than the systematic er- 
rors in the measurement. For more peripheral collisions, Raa^ increases with 
iVpart consistent with the modest suppression of high-pr production shown for 
peripheral collisions in the top panel of Fig. 38. The initial rise and subsequent 
decrease of Raa^ with increasing iV par t suggests that the high-p r hadron yield 
in Au+Au collisions has no simple dependence on iV pa rt- The observation that 
the high-pr yields initially increase proportional to Tab demonstrates that in 
the most peripheral Au+Au collisions the hard-scattering yields are consistent 
with point-like scaling. However, the deviation from Tab scaling sets in rapidly, 
becoming significant by iV pa rt = 50. By iV par t = 100 the high-p T suppression 
is so strong that high-p^ yields grow even more slowly than proportional to 

iVpart • 



6.2 xt scaling in Au+Au collisions at RHIC 




Fig. 39. Power-law exponent n(xx) for 7r° and h spectra in central and peripheral 
Au+Au collisions at y / ijv r /v = 130 and 200 GeV [53]. 
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If the production of high-p<r particles in Au+Au collisions is the result of hard 
scattering according to pQCD, then xt scaling should work just as well in 
Au+Au collisions as in p + p collisions and should yield the same value of 
the exponent ti(xt, y/s)- The only assumption required is that the structure 
and fragmentation functions in Au+Au collisions should scale, in which case 
Eq. 23 still applies, albeit with a G(xt) appropriate for Au+Au. In Fig. 39, 
n(xT, ^/snn) in Au+Au is derived from Eq. 23, for peripheral and central 
collisions, by taking the ratio of Ed 3 a/dp 3 at a given xt for ^/snn = 130 and 
200 GeV, in each case. The 7r°'s exhibit x T scaling, with the same value of 
n = 6.3 as in p+p collisions, for both Au+Au peripheral and central collisions, 
while the non-identified charged hadrons xr-scale with n = 6.3 for peripheral 
collisions only. Notably, the h^ 1 in Au+Au central collisions exhibit a signifi- 
cantly larger value of n, indicating different physics, which will be discussed 
below. The xt scaling establishes that high-p^ tt° production in peripheral and 
central Au+Au collisions and h^ 1 production in peripheral Au+Au collisions 
follow pQCD as in p + p collisions, with parton distributions and fragmenta- 
tion functions that scale with xt, at least within the experimental sensitivity 
of the data. 



6.3 Two-hadron azimuthal- angle correlations 

We argued in Sec. 5 that the production of hadrons at high-p T results pre- 
dominantly from hard scattering followed by fragmentation of the outgoing 
parton(s). While this result is well established in p{p) + p collisions, it might 
not be true in Au+Au collisions when the yield of high-pr particles is modi- 
fied so dramatically compared to expectations. Since a hard-scattered parton 
fragments into multiple particles within a restricted angular region (i.e. a jet) 
a reasonable way to check the assumption that high-p^ hadron production 
in Au+Au collisions is due to hard scattering is to directly observe the an- 
gular correlations between hadrons in the jets. None of the experiments at 
RHIC are currently capable of reconstructing jets in the presence of the large 
soft background of a Au+Au collision. However, both STAR [199,200] and 
PHENIX [201,202] have directly observed the presence of jets by studying 
two-hadron azimuthal-angle correlations. Figure 40 shows preliminary distri- 
butions [201] of the relative azimuthal angle (A</>) between pairs of charged 
particles detected within the PHENIX acceptance in <i+Au collisions and pe- 
ripheral (60-90%) and central (0-10%) Au+Au collisions after the subtraction 
of combinatoric background. The pairs of particles are chosen such that one 
particle lies within a "trigger" p T range (2.5 < pr trig < 4 GeV/c) while the 
other "associated" particle falls within a lower pt window 1.0 < px < 2.5 
GeV/c. The distributions show the differential yield per A0 of associated par- 
ticles per detected trigger particle within the given pt ranges and within the 
rj acceptance of the PHENIX central arms (—0.35 < rj < 0.35). 
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Fig. 40. Differential yields per A(f> and per trigger particle of pairs of charged 
hadrons in d+Au, peripheral Au+Au and central Au+Au collisions at s/snn = 200 
GeV[201]. The pairs were selected with the higher-momentum "trigger" particle in 
the range 2.5 < pt < 4.0 GeV/c and the lower-momentum "associated" particle in 
the range 1.0 < px < 2.5 GeV/c. A constant background has been subtracted for 
all three distributions. 

The peaks observed at A0 = (near side) reflect the correlation between 
hadrons produced within the same jet while the broader peaks observed at 
A0 = 7r (away side) reflect the correlations between hadrons produced in 
one jet and hadrons produced in the "balance" jet. In the Au+Au cases, a 
cos2A0 modulation underlies the jet angular correlations due to the elliptic 
flow of particles in the combinatoric background and possibly also in part due 
to azimuthal anisotropies in the jets themselves (see below). Nonetheless, the 
cos2A0 contribution has little effect on the narrow same-jet (near-side) peak 
in the A0 distribution. 

We observe that the angular widths of the same-jet correlations are the same 
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within errors in all three data sets in spite of the factor of two larger yield 
of associated hadrons in central Au+Au collisions compared to <i+Au and 
peripheral Au+Au collisions. This result is demonstrated more quantitatively 
in Fig. 41 which shows the centrality dependence of the Gaussian widths of 
the same-jet peaks in the Au+Au A<fi compared to the jet widths extracted 
from d+Au collisions [201]. We see that the Au+Au two-hadron correlation 
functions show peaks with the same jet width as <i+Au collisions. Since this 
width is a unique characteristic of the parton fragmentation process, we con- 
clude that high-p T hadrons in Au+Au collisions result from hard scattering 
followed by jet fragmentation regardless of any medium modifications of the 
fragmentation multiplicity. 



6.4 High-pT suppression and energy loss 

The suppression of the production of high-p^ hadrons in heavy ion collisions at 
RHIC had been predicted long before RHIC started running [173,174,175,176,177,178,179,203]. 
It is now generally accepted that partons propagating in colored matter lose 
energy predominantly through medium-induced emission of gluon radiation 
[204,205]. An energetic parton scatters off color charges in the high-parton- 
density medium and radiates gluon bremsstrahlung. The reduction in the par- 
ton energy translates to a reduction in the average momentum of the fragmen- 
tation hadrons, which, in turn, produces a suppression in the yield of high-pr 
hadrons relative to the corresponding yield in p + p collisions. The power-law 
spectrum for px > 3 GeV/c implies that a modest reduction in fragmenting 
parton energy can produce a significant decrease in the yield of hadrons at a 
given pt- Thus, the suppression of the yield of high-p^ hadrons is generally 
believed to provide a direct experimental probe of the density of color charges 
in the medium through which the parton passes [206,188,194]. However, before 
proceeding to an interpretation of our results, we briefly discuss the theoret- 
ical understanding of the radiative energy loss mechanism and limitations in 
that understanding. 

The dominant role of radiative gluon emission was identified early on [175], 
but it took several years and much effort before rigorous calculations of the en- 
ergy loss taking into account Landau-Pomeranchuk-Migdal suppression [176] 
and the time evolution of the medium were available. Initial estimates of the 
radiative energy loss suggested an approximately constant AE/Ax [174,176], 
but later calculations [177,207,180,181] showed that the quantum interference 
can produce a loss of energy that grows faster than linearly with the propaga- 
tion path length, L, of the parton in the medium. However, this ideal growth 
of AE/Ax with increasing path length is never realized in heavy ion collisions 
due to the rapid decrease of the energy density and the corresponding color 
charge density with time [179,188,194,185]. Generally, all energy loss calcula- 
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Fig. 41. The azimuthal angle width of jets in 200 GeV Au+Au collisions extracted 
as the cr's of Gaussian fits to the 0° peak in the two-charged-hadron azimuthal- angle 
(A(p) correlation functions [201]. The correlation functions were formed from pairs 
with trigger hadron in the px range 2.5 < pr < 4.0 GeV/c and the associated 
hadron in the range 1.0 < px < 2.5 GeV/c. The dashed lines show the ±1<7 range 
of the jet widths in <i+Au collisions using the same momentum bins. In the Au+Au 
data, the effect of the elliptic flow has been subtracted in the extraction of the jet 
width. 

tions predict that the fractional energy loss of a propagating parton decreases 
with increasing parton energy. However, the precise evolution with parton en- 
ergy depends on the assumptions in the energy loss models and on the treat- 
ment of details like kinematic limits and non-leading terms in the radiation 
spectrum [205,204]. There are many different calculations of medium-induced 
energy loss currently available based on a variety of assumptions about the 
thickness of the medium, the energy of the radiating parton, and the coher- 
ence in the radiation process itself (see [208,205,204] for recent reviews). The 
Pt dependence of the PHENIX n° Raa values has ruled out the possibility of 
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a constant (energy independent) AE/Ax [203] and the original BDMS energy 
loss formulation (which the authors argued should not be applied at RHIC en- 
ergies). In fact, the only detailed energy loss model that predicted the flat px 
dependence of Raa over the pt range covered by RHIC data was the GLV pre- 
scription [181,209,210,211,194]. In the GLV formulation, the fractional energy 
loss for large jet energies varies approximately as log(E)/E but the authors 
observe that below 20 GeV the full numerical calculation of the energy loss 
produces a nearly constant AE/E [205]. However, the same authors argue 
that the flat Raa(pt) observed at high p T at 200 GeV also requires an ac- 
cidental cancellation of several different contributions including the separate 
Pt dependences of the quark and gluon jet contributions, the pt dependence 
of the Cronin enhancement, and shadowing/EMC effect. A comparison of the 
GLV results for the px dependence of the 7r° suppression to the PHENIX data 
is shown in Fig. 42. 
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Fig. 42. Comparisons of energy loss calculations [212,185] used to extract estimates 
for the initial parton number or energy density (see text for details) to the central 
200 GeV Au+Au vr° R A a(pt) measured by PHENIX. The Wang curves compare 
results with and without energy absorption from the medium. 



One of the most critical issues in the energy loss calculation is the treatment 
of the time evolution of the energy density of the matter through which the 
radiating parton is propagating. Even if transverse expansion of the created 
matter is ignored, the longitudinal expansion produces a rapid reduction in the 
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energy density as a function of time. Most energy loss calculations assume that 
the color charge density decreases as a function of proper time as p(r) = Po 7 o/' r 
in which case the measured Raa can be used to infer the product po r o- Here 
t represents the formation time of the partons from which the medium is 
composed and po the initial number density of those partons. Since the gluons 
have the largest cross section for scattering with other partons, the initial 
color-charge density is interpreted as the gluon density. Making the usual 
assumption that the produced partons are spread over a longitudinal spatial 
width S z = r 5y, the GLV authors relate the product p r o to the initial dn g /dy 
and obtain dnjdy = 1000 ± 200 from the PHENIX vr° R AA values [194]. The 
sensitivity of the GLV calculations to the details of the description of the 
transverse parton density and the transverse expansion of the matter has been 
tested by using the results of hydrodynamic calculations of the energy density 
as a function of position and time [213]. The average energy loss for partons in 
central Au+Au collisions evaluated under dramatically different assumptions 
was shown to be remarkably insensitive to details of the description of the 
parton density. The GLV results are also potentially sensitive to a "screening 
mass" that determines both the transverse momentum distributions of the 
virtual gluons absorbed from the medium in the bremsstrahlung process and 
an energy cutoff for the radiated gluons. This mass is related to the local 
energy density using lattice QCD calculations of the plasma screening mass 
[194]. However, it was shown by the authors that a factor of two change in the 
screening mass produces only a 15% change in the dn g /dy needed to describe 
the data. 

An alternative analysis of parton energy loss [214] starts from explicit cal- 
culation of higher-twist matrix elements for e + A collisions that account for 
coherent rescattering of the struck quark in the nucleus. The contributions of 
these higher-twist terms can be incorporated into modified jet fragmentation 
functions, producing an effective energy loss. This calculation can reproduce 
[188] the HERMES measurements of modified jet fragmentation in nuclear 
deep-inelastic scattering [215]. By relating the modified fragmentation func- 
tions from the higher-twist calculation to energy-loss results obtained from 
the leading term in an opacity expansion calculation (e.g. GLV) of medium- 
induced energy loss the parameters describing the rescattering in the nucleus 
in e + A collisions can be related to the parameters describing the medium in 
an explicit energy-loss calculation. By relating the two sets of parameters, the 
parton density in the hot medium can be related to the parton density in a 
cold nucleus [188]. Results of this analysis are shown in Fig. 42 for parameters 
that give an initial energy loss per unit length of 13.8 ±3.9 GeV/fm when the 
HIJING [79] parameterization of shadowing is used [185] (Note: this result is a 
factor of two larger that in [188] which was based on analysis of the 130 GeV 
results). However, an alternative (EKS) [216] shadowing description results 
in an initial energy loss of 16.1 ± 3.9 GeV/fm [185] in the same calculation 
indicating at least a 25% systematic error in the energy loss estimates due 
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to uncertainties in the description of nuclear shadowing. Nonetheless, these 
initial-energy-loss values are much larger than the time-averaged energy loss 
extracted from the calculation, 0.85 ± 0.24 GeV/fm for HIJING shadowing 
[185], due to the assumed 1/r decrease in the color-charged density. In fact, 
the average energy loss per unit path length in central Au+Au collisions [188] 
is comparable to the value for cold nuclear matter extracted from HERMES 
data [188]. However, the initial energy loss is estimated by Wang to be a fac- 
tor of ~ 30 larger than that in a cold nucleus [185] implying that the initial 
Au+Au parton density is larger by a factor > 30 than in cold nuclear matter 
[217]. 

As shown in Fig. 42 the Wang higher-twist calculation predicts a suppression 
that varies strongly with p T over the range where the experimental Raa(pt) 
values are flat. However, Wang and Wang have argued that absorption of 
energy from the medium needs to be accounted for in calculating the energy 
loss of moderate-pT partons [187]. They provide a formula which incorporates 
both parton energy loss and "feedback" from the medium that can reproduce 
the shape of the observed high-p T suppression as shown by the lower curve 
in Fig. 42. This formula, then, provides the energy loss estimate given above. 
This explanation for the observed pr independence of Raa, a crucial feature 
of the experimental data, is disquieting, however, because it contradicts the 
explanation provided by the GLV model which provides a consistent estimate 
of the initial energy density. The feedback of energy from the medium is not 
included in the GLV calculations and if this contribution is significant, then 
the agreement of the GLV predictions with the 7r° Raa(pt) over the entire 
Pt range would have to be considered "accidental". Also, the variation of the 
suppression in the Wang higher-twist calculation with p? reflects the AE oc 
log E variation of parton energy loss naturally obtained from approximations 
to the full opacity expansion [205]. As noted above, the GLV approach finds 
that incorporating non-leading terms in the opacity expansion produces AE oc 
E. Thus, while the absorption of energy from the medium in the Wang et al. 
approach may only be significant below pt = 5 GeV/c, the differences between 
the variation of energy loss with parton energy in the two approaches will not 
be confined to low p T . 

One source of uncertainty in the interpretation of the high-p-^ suppression is 
the role of possible inelastic scattering of hadrons after fragmentation. It was 
originally argued that final-state inelastic scattering of hadrons could produce 
all of the observed suppression [218]. The persistence of the jet signal with 
the correct width in Au+Au collisions would be difficult to reconcile with this 
hypothesis. Indeed, more recent analyses [219] discount the possibility that 
hadronic re-interaction could account for the observed high-pr suppression 
and indicate that only ~ 1/3 of fragmentation hadrons undergo final-state 
inelastic scattering [219]. Wang has also argued [220] that the complete pat- 
tern of high-p T phenomena observed in the RHIC data cannot be explained 
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by hadronic rescattering. However, this leaves open the question of whether 
hadronic re-interactions after jet fragmentation can be partially responsible 
for the observed high-pp suppression. There are a number of other open is- 
sues with the quantitative interpretation of the observed high-p T suppression. 
The calculations all assume that the jets radiate by scattering off static color 
charges while the typical initial gluon is often assumed to be ~ 1 GeV. 
Also the radiated gluons are assumed to be massless though a plasmon cutoff 
equal to the screening mass is applied. The systematic errors introduced by 
these and other assumptions made in the current energy loss calculations have 
not yet been evaluated though the gluon screening mass is being included in 
analyses of heavy-quark energy loss. 



6. 5 Empirical Energy Loss Estimate 



The observation that the suppression of high-p^ particle production is ap- 
proximately independent of pt above 4 GeV/c and that the p + p pr spectra 
are well described by a pure power-law function in the same px range allows 
a simple empirical estimate of the energy loss of hard-scattered partons in 
the medium. The 7r° invariant cross section measured by PHENIX in p + p 
collisions [60] is found to be well described by a power law 

d 3 n 1 d 2 n A 

rp f98^ 

dp 3 2ir prdprdy pr n 



for pr > 3.0 GeV/c with an exponent n — 8.1 ± 0.1. If we assume that none 
of the hard-scattered partons escape from the medium without losing energy, 
then the approximately p^-independent suppression above 4.5 GeV/c can be 
interpreted as resulting from an average fractional shift in the momentum of 
the final-state hadrons due to energy loss of the parent parton. The suppressed 
spectrum can be evaluated from the unsuppressed (p + p) spectrum by noting 
that hadrons produced in Au+Au collisions at a particular px value, would 
have been produced at a larger pt value pt = Pt + S(pt) in p + p collisions. 
If the energy loss is proportional to pr then we can write S(pr) = So£>t so 
Pt' = (1 + S )pt Then, the number of particles observed after suppression in 
a given Ap T interval is given by 

dn dn dpr' A 
~d~p~ T = dp T ' dp T = (1 + So)^ Pt^ ' 



We note that the factor accounts for the larger relative density of parti- 
cles per measured px interval due to the effective compression of the pr scale 
caused by the induced energy loss; this factor is necessary for the total num- 
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ber of particles to be conserved. The nuclear modification factor then can be 
expressed in terms of So, 

= (i + s^i - (30) 



Using this very simple picture, we can estimate the fraction of energy lost by 
hard-scattered partons in the medium from our measured R AA values. First 
we obtain S from Eq. 30 

So = rJiw " L (31) 



Then we observe that the hadrons that would have been produced in p + p 
collisions at a momentum (1 + So)pr were actually produced at pr, implying 
a fractional energy loss 

S, oss = 1 - 1/(1 + S ) = 1- R A A l/{n - 2) - (32) 



Figure 43 shows the centrality dependence of S\ oss obtained from the px- 
averaged Raa values shown in Fig. 38. For the most central Au+Au collisions 
at 200 GeV we obtain S\ oss = 0.2, which naively implies that an average 20% 
reduction in the energy of partons in the medium will produce the suppres- 
sion observed in the 7r° spectra above 4.5 GeV/c. The extracted S\ oss values 
are well described by an iVp art 2 / 3 dependence using the most central bin to fix 
the proportionality constant. This result agrees with the GLV prediction for 
the centrality dependence of the medium-induced energy loss. 

It has been shown previously [221,222] that fluctuations in the radiation pro- 
cess can distort an estimate of parton energy loss using the procedure de- 
scribed above. Because of the steeply falling px spectrum, the partons that 
lose less energy dominate the yield at a given px so our determination of S\ oss 
will significantly underestimate the true energy loss. However, it has also been 
observed that this distortion can largely be compensated by a single multi- 
plicative factor of value ~ 1.5 — 2 [221]. While we cannot use the empirically 
extracted energy loss to estimate an initial gluon density, we can evaluate 
the consistency of our results with estimates of (dE/dx) in the medium. If 
we take into account the factor of 1.5 — 2.0 renormalization of the S\ oss , we 
estimate that 10 GeV partons lose ~ 3 — 4 GeV of energy. If the typical path 
length of these partons is on the order of the nuclear radius then we can in- 
fer a AE/Ax ~ 0.5 GeV/fm which is in good agreement with the estimate 
from Wang [188]. We can also use the above empirical energy loss approach to 
evaluate possible systematic errors in the estimate of the initial gluon density. 
For example, if one third of the observed suppression were a result of final- 
state hadronic interactions in the medium, then the suppression due to energy 
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Fig. 43. Calculated energy loss shift factor, S\ oss vs. N pa _ Tt for ir° and charged hadron 
production in 200 GeV Au+Au collisions. The band around the values indicates 
systematic errors resulting from uncertainties in Tab an d the normalization of the 
p + p spectrum. The dot-dashed curve shows an iV par t 2 ' 3 scaling of S\ oss using the 
most central bin to fix the proportionality constant. 



loss would be a factor of 1.5 smaller than that implied by the measured Raa 
values, assuming that every fragmentation hadron that interacts effectively 
"disappears" by being shifted to much lower momentum. As a result, Sq in 
central Au+Au collisions would be reduced from 0.25 to 0.17, implying 30% 
reduction in the estimated energy loss. If the energy loss is indeed proportional 
to the initial gluon density then the uncertainty in the effect of the final-state 
hadronic interactions would introduce a 30% systematic error in dn g /dy. 
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6.6 Conclusions 



The observed suppression of high-p T particle production at RHIC is a unique 
phenomenon that has not been previously observed in any hadronic or heavy 
ion collisions at any energy. The suppression provides direct evidence that 
Au+Au collisions at RHIC have produced matter at extreme densities, greater 
than ten times the energy density of normal nuclear matter and the highest 
energy densities ever achieved in the laboratory. Medium-induced energy loss, 
predominantly via gluon bremsstrahlung emission, is the only currently known 
physical mechanism that can fully explain the magnitude and pt dependence 
of the observed high-p^ suppression. This conclusion is based on evidence 
provided above that we summarize here: 

• Observation of the xt scaling of the high-p T hadron spectra and measure- 
ments of two-hadron azimuthal-angle correlations at high pt confirm the 
dominant role of hard scattering and subsequent jet fragmentation in the 
production of high-p T hadrons. 

• d+Au measurements demonstrate that any initial-state modification of nuclear- 
parton distributions has little effect on the production of hadrons with 
Pt > 2 GeV/c at mid-rapidity. 

• This conclusion is further strengthened by preliminary PHENIX measure- 
ments showing that the yield of direct photons with pt > 5 GeV/c is consis- 
tent with a Tab scaling of a pQCD-calculated p + p direct-photon spectrum. 

• Analyses described above indicate that final-state hadronic interactions can 
only account for a small fraction of the observed high-px suppression. 

Interpreted in the context of in-medium energy loss, the high-pr suppression 
data rule out the simplest energy loss prescription — a jet energy indepen- 
dent AE/Ax. The approximately flat Raa(pt) was predicted by the GLV 
energy loss model from which the most explicit estimates of the initial gluon- 
number density, dn g /dy = 1000 ± 200 and a corresponding initial energy den- 
sity Eq ~ 15GeV/fm 3 [194], have been obtained. An alternative estimate from 
the analysis of Wang et al. [188] yields a path-length-averaged energy loss of 
0.5 GeV/fm. Assuming a 1/r time evolution of the energy density a much 
larger initial energy loss of 13-16 GeV/fm is obtained. That estimate com- 
bined with the estimated 0.5 GeV/fm energy loss of partons in cold nuclear 
matter yields an initial Au+Au gluon density > 30 times larger than that in 
nuclei [217]. From this result, Wang concludes that the initial energy density is 
a factor of ~ 100 times larger than that of a nucleus which would correspond 
to 16GeV/fm 3 [217]. While this conclusion is consistent with the independent 
estimate from GLV, we note that the two models provide completely different 
explanations for the nearly p^-independent Raa - the most unique feature 
of the single-particle high-pp suppression - and the differences between the 
approaches may not be confined to low p T . An empirical analysis of the par- 
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ton energy loss suggests that the Wang estimate of > 0.5 GeV/fm for the 
average parton AE/Ax is consistent with the measured Raa values in central 
Au+Au collisions. However, some outstanding issues with current energy loss 
calculations and the interpretation of high-p T suppression were noted above. 
Most notably, rescattering of hadrons after parton fragmentation could affect 
the observed high-pr suppression even if such rescattering cannot explain the 
pattern of jet quenching observations. Using results from [219] and our empir- 
ical energy loss analysis, we estimated that hadronic interactions could modify 
extracted values for initial parton densities by only 30%. However, we cannot 
evaluate the potential systematic error in extracted parton densities due to 
other untested assumptions of the energy loss calculations. Therefore, to be 
conservative we interpret the extracted initial gluon number and energy den- 
sities as order-of-magnitude estimates. Even then, the 15GeV/fm 3 estimated 
by Gyulassy and Vitev from the central 200 GeV Au+Au ir° Raa(pt) mea- 
surements indicates that the matter produced in central Au+Au collision has 
an energy density > 10 times normal nuclear matter density. 



7 HADRON PRODUCTION 



Descriptions of heavy ion collisions have provided an understanding of early 
energy densities, production rates and medium effects of hard partons, and col- 
lective flow of matter. However, hadronization — the process by which partons 
are converted into hadrons — is not well understood. The process of hadroniza- 
tion is particularly important since it includes both the dressing of the quarks 
from their bare masses, i.e. the breaking of approximate chiral symmetry, 
and the confinement of quarks into colorless hadrons. One could conclude 
that a quark-gluon plasma had been formed if one had conclusive evidence of 
hadronization occurring from a thermal distribution of quarks and gluons. 

Hadronization processes have been studied over many years in proton-proton 
and electron-positron reactions. Hadron formation, by its very nature a non- 
perturbative process, has often been parameterized from data (e.g. fragmen- 
tation functions D(z)) or phenomenologically described (e.g. string models) 
[223]. From QCD one expects that hadron production at high transverse mo- 
mentum is dominated by hard scattering of partons followed by fragmentation 
into "jets" or "mini-jets" of hadrons. Following the assumptions of collinear 
factorization, the fragmentation functions should be universal. This univer- 
sality has proved a powerful tool in comparing e + e~ annihilation to hadron- 
hadron reactions. One feature of jet fragmentation is that baryons and an- 
tibaryons are always suppressed relative to mesons at a given px [224,225]. 
Phenomenologically this can be thought of as a large penalty for creating a 
diquark-antidiquark pair for baryon formation vs. a quark-antiquark pair for 
meson formation. 
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In hadron-hadron reactions, hard scattering followed by fragmentation is con- 
sidered to be the dominant process of hadron production for particles with 
Pt > 2 GeV/c at mid-rapidity. At low transverse momentum, where particles 
have p T < 2 GeV/c, particle interactions are often referred to as "soft". In 
small momentum transfer reactions the effective wavelength of interactions is 
longer than the spacing of individual partons in a nucleon or nucleus. Thus 
coherence effects are expected to result in large violations of factorization and 
universality of fragmentation functions. Hadron formation mechanisms in this 
"soft" regime are poorly understood. We are particularly interested in the 
study of hadron formation in the region of p T pa 2-5 GeV/c, where production 
is expected to make the transition from "soft" to "hard" mechanisms. 

7. 1 Baryons and Antibaryons 

One of the most striking and unexpected observations in heavy ion reactions 
at RHIC is the large enhancement of baryons and antibaryons relative to pions 
at intermediate p T pa 2-5 GeV/c. As shown in Fig. 44, the (anti)proton to pion 
ratio is enhanced by almost a factor of three when one compares peripheral 
reactions to the most central gold-gold reactions [52] 7 . This of course is in 
sharp contrast to the suppression of pions in this region. 

We can investigate this (anti)baryon excess to much higher pr by comparing 
our inclusive charged spectra (primarily pions, kaons and protons) with our 
neutral pion measurements [52]. Shown in Fig. 45 is the charged hadron to 
7r° ratio as a function of transverse momentum in ten centrality bins. We 
observe a significant increase of the {h + + h~)/ir° ratio above 1.6 in the pr 
range 1-5 GeV/c that increases as a function of collision centrality. The ratio 
of h/n — 1.6 is the value measured in p + p reactions [224], and is thought 
to arise from jet fragmentation. In Au+Au central reactions, above pt ~ 5 
GeV/c, h/ir returns to the p + p measured baseline. This implies that the 
(anti)baryon excess occurs only in the limited pr window pa2-5 GeV/c, and 
then returns to the universal fragmentation function expectation. 

As discussed in section 6, pions in this pt range are suppressed by almost a 
factor of five relative to binary collision scaling for central Au+Au reactions. 
Thus, one possible interpretation of the large (anti)proton to pion ratio is that 
somehow the baryons are not suppressed in a manner similar to the pions. 
Figure 46 shows that in fact (anti)proton production appears to follow binary 
collision scaling over the transverse momentum range pt = 2-5 GeV/c [52]. 
However, the h/n° ratios shown in Fig. 45 imply that above pr > 5 GeV/c, 
the (anti)protons must be as suppressed as the pions. 

7 All PHENIX (anti)proton spectra shown in this section are corrected for feed 
down from heavier resonances. 
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Fig. 44. p/ir (left) andp/ir (right) ratios for central (0-10%), mid-central (20-30%) 
and peripheral (60-92%) Au+Au collisions at ^/sWn = 200 GeV[52]. Open (filled) 
points are for 7r + /~ (t^ ), respectively. Data from ^/s = 53 GeV p+p collisions [224] 
are shown with stars. The dashed and dotted lines are (p + p)/(tt + + 7t~) ratio in 
gluon and in quark jets [225]. 



Characteristics of the intermediate pt (anti)protons are: 



• A large enhancement of the p/n and p/ir ratios in central Au+Au collisions. 

• A ratio in peripheral collisions which is in agreement with that from p + p 
collisions. 

• A smooth increase from peripheral to central Au+Au collisions. 

• A similar effect for protons and antiprotons. 

• Approximate scaling of (anti)proton production at px ~ 2-4 GeV/c with 
the number of binary nucleon-nucieon collisions. 

• Suppression relative to binary collision scaling similar for (anti)protons and 
pions for px > 5 GeV/c. 

Large proton to pion ratios have also been observed in heavy ion collisions at 
lower energies. Figure 47 shows pt distributions of protons, antiprotons, and 
pions in central Pb+Pb collisions at the SPS and in central Au+Au collisions 
at the AGS. The p/n ratio in central Pb+Pb collisions at the SPS is greater 
than unity for p T > 1.3 GeV/c. At the AGS, the proton spectrum crosses 
pion spectra at px ~ 0.5 GeV/c, and the p/n ratio is about 20 at pt = 1-6 
GeV/c. The p/V ratios in the low-energy heavy- ion collisions are also enhanced 
compared with p + p collisions at the same energy. 
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Fig. 45. Charged hadron to tt° ratio for different centrality classes for Au+Au colli- 
sions at y/SNN = 200 GeV[53]. Error bars represent the quadratic sum of statistical 
and point to point systematic errors. The shaded band shows the normalization 
error common to all centrality classes. The line at 1.6 is the h/ir ratio measured in 
p + p collisions at ^fs = 53 GeV [224] and e+e- collisions [225] . 
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Fig. 46. p and p invariant yields scaled by N co u in Au+Au collision at ^snn = 200 
GeV[52]. Error bars are statistical. Systematic errors on N co u range from ~10% for 
central to ~28% for 60-92% centrality. Multiplicity dependent normalization errors 
are «3%. 



Most of the protons in these lower-energy heavy-ion collisions are not produced 
in the collision. Rather they are protons from the beam or target nucleus (Pb 
or Au) that are transported to large px at mid-rapidity. As discussed in section 
3, a strong radial flow with velocity (3t ~ 0.5 is produced in heavy ion collisions 
at AGS and SPS energies. The large p/n ratio can be interpreted as a result of 
this radial flow. Since the proton is heavier, a fixed velocity boost results is a 
larger momentum boost than for pions, and thus enhances p/w ratio at higher 
Pt- In contrast, at RHIC energies, most of protons are produced particles [40]. 
The anomalously large antibaryon-to-meson ratio p/ir ~ 1 at high pt > 2 
GeV/c is a unique result from RHIC. Such a large p/n ratio has not been 
observed in any other collision system. Figure 47 shows that p/n is less than 
~ 0.1 at the SPS, and it is less than 1/100 at the AGS. It should also be noted 
that the measurements from the AGS/SPS are limited to lower p? {pt < 2 
GeV/c), where soft physics is still dominant, while at RHIC we observe a large 
p(P)/tt ratio in p T f& 2-5 GeV/c where hard processes are expected to be the 
dominant mechanism of particle production. 
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Fig. 47. Invariant yields of p, p, and tt as function of px in central Pb+Pb collisions 
at the SPS (^/snn = 17 GeV) (left panel) and in central Au+Au collisions at the 
AGS (y / ijv]v = 5 GeV) (right panel). The p spectrum from the AGS is scaled up 
by a factor 100. All data are at mid-rapidity {y — y C m ~ 0) and are from W98 [163], 
NA44 [226], NA49 [227], and E866 [228,229]. 



7.2 The (f) Meson 



We have extended our identified hadron studies to include the <p vector meson 
as measured in the K + K~ decay channel. The is a meson, and is in that 
sense similar to the pion with a valence quark and antiquark, and yet its mass 
is comparable to that of the proton. 

Figure 48 shows Rcp, the ratio of production in central to peripheral Au+Au 
collisions scaled by binary collisions, for protons, pions and <ft mesons detected 
via its KK decay channel [69] in Au+Au collisions at ^snn = 200 GeV. A 
large suppression of pions at pr > 2 GeV/c is observed (as detailed in Section 
6), and a lack of suppression for the protons and antiprotons as expected 
from Fig. 46. The <ft follows the suppression pattern of the pions within errors, 
indicating that the surprising behavior of the protons is not followed by the <fi. 
Figure 49 shows a comparison between the pt spectral shape for protons and 
the in central and peripheral Au+Au reactions. The two spectra agree with 
each other within errors for the most central events. Thus, although the yields 
are evolving differently with collision centrality, giving rise to the deviation 
from unity of Rcp, the pr distributions appear quite similar. 
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Fig. 48. The Rcp of the <p as measured in the KK channel, compared to the protons 
and pions for Au+Au collisions at ^Jsnn = 200 GeV[69]. 

7.3 Jet Correlations 

A crucial test of the origin for the enhanced (anti)proton to pion ratio is to 
see if baryons in this intermediate pt regime exhibit correlations characteristic 
of the structure of jets from hard-scattered partons. Particles which exhibit 
these correlations are termed "jet-like". Figure 50 shows the associated partner 
particle yield within the relative angular range 0.0 < <p < 0.94 radians on the 
same side as trigger baryons and mesons [71] . Correlated pairs are then formed 
between the trigger particle and other particles within the above mentioned 
angular range. Mixed events are used to determine the combinatorial (i.e. 
non-jet-like) background distribution, which is subtracted after modulation 
according to the measured v<i- 

The partner yield increases for both trigger baryons and mesons by almost 
a factor of two from deuteron-gold to peripheral and mid-central Au+Au re- 
actions. We then observe a decrease in the jet-like correlations for baryons 
relative to mesons for the most central collisions. It is notable that this obser- 
vation is of limited significance within our current statistical and systematic 
errors. Over a broad range of centrality 10-60% the partner yield is the same 
for protons and pions within errors. This is notable since the (anti)proton to 
pion ratio has already increased by a factor of two for mid-central Au+Au 
relative to proton-proton reactions, with the implication that the increase in 
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Fig. 49. (p+p) and invariant yield as a function of transverse momentum for central 
0-10% and peripheral 40-92% Au+Au reactions at ^/snn = 200 GeV. The two 
distributions for each centrality class are given an arbitrary relative normalization 
to allow for comparison of the pt dependent shapes. Data are from [69]. 

the p/n ratio is inclusive of the particles with jet-like correlations. 

The dashed line in Figure 50 shows the expected centrality dependence of 
partners per baryon if all the "extra" baryons which increase the p/n over 
that in p+p collisions were to arise solely from soft processes. Baryons from 
thermal quark recombination should have no jet-like partner hadrons and 
would dilute the per-trigger conditional yield. Because this simple estimate 
does not allow for meson production by recombination, which must also occur 
along with baryon production, it represents an upper limit to the centrality 
dependence of jet partner yield from thermal recombination. The data clearly 
disagree with both the centrality dependence and also the absolute yields of 
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this estimation, indicating that the baryon excess has the same jet-like origin 
as the mesons, except perhaps in the highest centrality bin. 
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Fig. 50. Centrality dependence of associated charged hadron yield (1.7 < px < 2.5 
GeV/c) above combinatorial background for trigger baryons and trigger mesons in 
the pt range 2.5-4.0 GeV/c in a 54° cone around the trigger particle in Au+Au 
collisions at -y/J/vTV = 200 GeV[71]. The error bars are statistical errors and the 
gray boxes are systematic errors. The dashed line represents an upper limit of the 
centrality dependence of the near-side partner yield from thermal recombination 
(see text). 

The characteristics of the jet-like particles are compared to inclusive hadrons 
in Fig. 51, which shows the centrality dependence of the pt distributions of 
jet-like partners and inclusive hadrons. One can see that, within the statistics 
available, the slopes of the associated particle spectra in p + p, d+Au, periph- 
eral and mid-central Au+Au collisions are very similar for both trigger mesons 
and trigger baryons. The partner spectra are harder than the inclusive hadron 
spectra, as expected from jet fragmentation. In the most central collisions, the 
number of particles associated with trigger baryons is very small, resulting in 
large statistical error bars. However, the inverse slopes of the jet-like partners 
and inclusive hadron distributions agree better in central collisions than in 
peripheral collisions. 

We can then make the following general observations: 

• Trigger (anti)protons and mesons have comparable near-side associated- 
particle yields over a broad range in centrality, indicating a significant jet- 
like component for both. 
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• There is an indication that the proton partner yield tends to diminish for 
the most central collisions, unlike for leading mesons. 

• Within the limited statistics available for the measurement, the inverse 
slopes of the associated particles are similar for both mesons and baryons. 
These are harder than for the inclusive spectra. 

• Trigger particles in Au+Au collisions appear to have more associated par- 
ticles than in d+Au collisions. This is true for all centralities aside from the 
most peripheral, and except for leading baryons in central collisions. 
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Fig. 51. The inverse slopes for the momentum distributions of the associated parti- 
cles shown in Fig. 50. The gray band is the inverse slope of the momentum distri- 
bution of the inclusive hadrons from Au+Au collisions at y^STVTV = 200 GeV[71]. 



7.4 Soft Physics 



In hadron-hadron reactions, hard scattering followed by fragmentation is con- 
sidered to be the dominant process of hadron production with pt > 2 GeV/c 
at mid-rapidity. However, as detailed in section 3, there is strong evidence 
for explosive collective motion of particles in the medium. If the mean free 
path for particles in the medium is small, then all particles must move with 
a common local velocity as described by hydrodynamics. Therefore, heavier 
particles receive a larger momentum boost than lighter particles. This effec- 
tive shifting of particles to higher p^ results in a "shoulder-arm" shape for the 
(anti)proton p T spectra, visible in Fig. 49. 
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7.4-1 Hydrodynamics 



Is it possible that this soft hadron production extends to higher p? for baryons 
than mesons? Hydrodynamic boosting of "soft" physics for heavier particles 
into the p T > 2 GeV/c offers a natural explanation for the enhanced p/n and 
p/n ratios [106]. 

As seen in Section 3, some hydrodynamical models can describe both the pro- 
ton and the pion spectra. Consequently, the p/n ratio is also reproduced (Fig. 
52). It is clear that the description of the p/n ratio is not unique and differ- 
ent calculations yield quite different results. Above some p?, hydrodynamics 
should fail to describe the data and fragmentation should dominate. Pure hy- 
drodynamics predicts that this ratio would continue to increase essentially up 
to pt — > oo. However, these particles cannot have a zero mean free path in the 
medium. Any finite mean free path and a finite volume will limit the number 
of pt "kicks" a particle can receive. For this reason many of the hydrodynamic 
calculations are not extended into the px region 2-5 GeV/c in which we are 
interested. 

Hydrodynamic calculations do not specify the quanta that flow; rather they 
assume an equation of state. When applied at RHIC, most calculations start 
with a quark-gluon-plasma equation of state and transition to a resonance gas. 
The mapping of the fluid onto hadrons is somewhat ad hoc, and often uses 
the Cooper- Frye freezeout [126], giving the typical hierarchies of momenta one 
sees where heavier particles receive a larger boost. 

As mentioned previously, this generic feature of a transverse velocity boost 
yielding an increase in the baryon to meson ratio relative to proton-proton 
reactions is not unique to RHIC as shown in Fig. 47. However, a major differ- 
ence between lower-energy results and those at v /saT/v = 200 GeV is that at 
these highest energies there is a significant hard-process contribution. If the 
source of the excess baryons is the transport of soft baryons to the intermedi- 
ate pt range, then it is purely coincidental that the baryons scale with binary 
collisions. More importantly, we should expect a significant decrease in the 
jet-like partner yield for baryons relative to mesons. Although there may be a 
hint of this for the most central reactions, one expects this decrease to follow 
the centrality dependence of the increase in p/n ratio. Thus, this effect should 
already reduce the partner yield by a factor of two in mid-central Au+Au 
reactions. This is ruled out by the data. 

7-4-2 Recombination Models 

The quark recombination or coalescence model is a different physics frame- 
work in which baryons receive a larger pt boost than mesons. These models 
were frequently invoked in the 1970's [230,231] in an attempt to describe the 
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Fig. 52. p/n ratios for central (0-10%) Au+Au collisions at ^/snn = 200 GeV[52] 
compared to hydrodynamic models [106,107,117,118,119]. 

rapidity distribution of various hadronic species in hadron — hadron reac- 
tions. More recently, these models have been applied to describe the forward 
charm hadron production in hadron — nucleus reactions at Fermilab [232]. 
In this case they calculate a significant probability for D meson formation 
from a hard-scattering-created charm quark with a light valence quark in the 
projectile. The quark coalescence mechanisms have some similarities to light 
nuclei coalescence. However, wave functions are relatively well determined for 
light nuclei, whereas the hadron wave functions are neither easily described 
by partons nor directly calculable from QCD. 



Recently, quark recombination has been successfully applied to describe a 
number of features of heavy ion collisions [233,234] (Duke model). In this 
picture, quarks in a densely populated phase space combine to form the final- 
state hadrons. This model uses the simplifying assumption that the mass is 
small relative to the momentum giving a prediction largely independent of the 
final hadron wave function 8 . The coalescing parton distribution was assumed 
to be exponential, i.e. thermal, and recombination applied for hadrons where 
m 2 /pT << 1. At very high particles are assumed to arise from fragmen- 
tation of hard partons with a standard power law distribution; the relative 
normalization of the thermal source with respect to this process is an impor- 
tant external parameter to the model. A crucial component of recombination 

8 The recombination model prediction of these models is independent of the final 
hadron wave function with an accuracy of about 20% for protons and 10% for pions 
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models is the assumption that the partons which recombine carry a mass 
which is essentially equal to the mass of the dressed constituent quarks 9 . If 
all observables of intermediate pt hadrons can be explained by recombina- 
tion of only thermal quarks, this would essentially prove the existence of a 
quark-gluon plasma in the early stage of the collisions. 

Three essential features are predicted by recombination models. First, baryons 
at moderate pr are greatly enhanced relative to mesons as their transverse mo- 
mentum is the sum of 3 quarks rather than 2. Recombination dominates over 
parton fragmentation in this region, because, for an exponential spectrum re- 
combination is a more efficient means of producing particles at a particular 
Pt- This enhancement should return to its fragmentation values at higher p?- 
In the intermediate range, all mesons should behave in a similar manner re- 
gardless of mass, as should all baryons. Secondly, recombination predicts that 
the collective flow of the final-state hadrons should follow the collective flow 
of their constituent quarks. Finally, recombination causes thermal features to 
extend to higher transverse momentum, p? » Tq than one might naively 
expect since the underlying thermal spectrum of the constituents gets a mul- 
tiplication factor of essentially 3 for baryons and 2 for mesons. A last general 
feature which is true for the simplest of the models, but may not necessarily be 
true for more complex models, is that at intermediate pr, recombination is the 
dominant mechanism for the production of hadrons — particularly of baryons. 

Other recombination calculations have relaxed the assumptions previously de- 
scribed, at the cost of much more dependence on the particular form of the 
hadronic wave function used. One such calculation [235,236] (Oregon model) 
uses a description of hadronization which assumes that all hadrons — including 
those from fragmentation — arise from recombination. Hard partons are al- 
lowed to fragment into a shower of partons, which can in turn recombine— 
both with other partons in the shower and partons in the thermal background. 
Another model [237] (TAMU model) uses a Monte-Carlo method to model the 
production of hadrons allowing recombination of hard partons with thermal 
partons, and includes particle decays, such as p — > 2n which produces low-]>r 
pions. 

Figure 53 shows several recombination model calculations compared to the pj tt 
ratio from PHENIX. The general features at pr > 3 GeV/c are reasonably 
reproduced — that is the protons show a strong enhancement at moderate Pt 
which disappears at pt > 5 GeV/c consistent with the measured h/n ratio 
shown in Fig. 45. The more complicated models do a better job, as one might 



9 The actual source of this mass is under discussion. It may be that the chiral 
phase transition is slightly above the deconfmement transition. In this case, the 
mass would be from the dressing of the quarks. Another possibility is that the mass 
is a thermal mass which happens to be similar to the constituent quark mass. 
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Fig. 53. The proton to pion ratio measured by PHENIX for Au+Au collisions at 
\J s nn = 200 GeV[52]. Several comparisons to recombination models as mentioned 
in the text are shown. 

expect in the px < 3 GeV/c region, where the assumptions made by the 
Duke model begin to break down. Since the recombination model's essential 
ingredient is the number of constituent quarks in a hadron, the similarity of 
Rcp for the (f> and pions is nicely explained. 

Figure 54 shows the fraction of hadrons arising from recombination of only 
thermal quarks, as a function of px- For px between 2.5 and 4 GeV/c the 
fraction of protons from recombination is greater than 90% for all impact pa- 
rameters, and is essentially 100% for the most central collisions. For pions the 
value is between 40 and 80%, depending on the centrality. This is contradicted 
by the data in Fig. 50 which clearly shows jet-like correlations for both pions 
and protons in mid-central collisions. It should be noted that the yield of par- 
ticles associated with baryons in very central collisions appears to decrease, 
indicating a possible condition where the simple picture of recombination of 
purely thermal quarks may apply. 

One can examine the general prediction for the elliptic flow of identified par- 
ticles by rescaling both the t>2 and the transverse momentum by the number 
of constituent quarks as shown in Fig. 55. This scaling was first suggested 
by Voloshin [238]. Above px/n of 1 GeV/c (corresponding to 3 GeV/c in the 
proton transverse momentum) all particles essentially plateau at a value of 
about 0.35 presumably reflecting the elliptic flow of the underlying partons. 
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Fig. 54. The ratio t{Pt) = R/ {R + F) of recombined hadrons to the sum of recom- 
bination (R) and fragmentation (F) for pions (solid), K°s (dashed) and p (dotted 
lines) [234] in Au+Au collisions at y / Sjv r /v = 200 GeV. For protons and pions differ- 
ent impact parameters 6 = 0, 7.5 and 12 fm (from top to bottom) are shown. K°s 
is for b = fm only. 

Interestingly, even at lower values of the transverse momentum, all particles 
also fall on the same curve aside from pions. 

It is clear from the jet correlations observed that the majority of moderate 
Pt baryons in peripheral and mid-central collisions cannot arise from a purely 
thermal source, as that would dilute the per-trigger partner yield. The jet 
structure and collision scaling indicate that at least some of the baryon excess 
is jet-like in origin. The relatively short formation time for baryons of such 
momenta suggests that allowing recombination of fragmentation partons with 
those from the medium may solve the problem and better reproduce the data. 
Both the Oregon and TAMU models have mechanisms to do this. However, 
such modification of the jet fragmentation function must also modify the ellip- 
tic flow, and could break the quark scaling needed to reproduce the observed 
t> 2 trends. Hence, the jet structure of hadrons at 2-5 GeV/c px presents a 
challenge to models of the hadron formation. 

Figure 56 shows a comparison of the elliptic flow calculated by the TAMU 
model [241] with PHENIX data from Au+Au collisions at ^/JJ^ = 200 GeV. 
The model includes the recombination of hard and soft partons, as well as the 



100 



c 0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



▲ 

★ 

□ 

o 

A 



PHENIX 7T+K" 20-40% 
STAR 7i, 11-45% 
STAR K°, 5-30% 
PHENIX ti° 20-40% 
PHENIX p, 20-40% 
STAR p, 11-45% 
STAR A, 5-30% 



f. 



A-L 



£IA 



^A 

A 



0.2 0.4 0.6 0.8 1 



1.2 1.4 1.6 1.8 2 

p T /n (Gev/c) 



Fig. 55. «2 as a function of transverse momentum for a variety of particles for Au + 
Au collisions where both v 2 and pr have been scaled by the number of constituent 
quarks in the particle. The meson data are shown with filled symbols; n~ + K~ 
from PHENIX at = 200 GeV [50] (filled circles), charged vr from STAR at 

y/I^N = 130 GeV [103] (filled squares), K° s from STAR at ^/IWn = 200 GeV [239] 
(filled triangles), and vr° from PHENIX at ^sWJj = 200 GeV [240] (filled stars). 
While the baryons are shown with open symbols; p from PHENIX at ^snn = 200 
GeV [50] (open squares), p from STAR at ^/snn = 130 GeV [103] (open circles), 
and A from STAR at = 200 GeV [239] (open triangles). 



decay of resonances such as the p. In this model, at least, the agreement of v 2 
with the data is preserved — in addition a simple explanation is given for the 
excess of pion v 2 at low pr- A similar conclusion in shown in [242]. This would 
seemingly attribute all the elliptic flow to the partonic phase leaving no room 
for additional flow to be produced in the later, hadronic stage — which may 
be in contradiction to hydrodynamic interpretations of the hadronic state as 
demanded by a variety of signatures such as the pt spectra of the protons 
and pions (see Section 3). It is clear that a more comprehensive comparison 
of observables should be undertaken to check the validity of these models. 
Higher-statistics jet studies with different identified particles by PHENIX in 
Run-4 will help clarify the situation. 
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Fig. 56. V2/n in the TAMU model, where n is the number of constituent quarks in 
a particle for protons and pions[241]. Scaled pion (dashed line) and proton (dotted 
line) results from the TAMU model are shown in addition to charged pion and 
proton measurements from the PHENIX experiment from minimum bias Au+Au 
reactions at ^snn = 200 GeV [50]. This model allows for the recombination of 
hard partons and soft partons, as well as the decay of resonances such as the two 
pion decay of the p meson. One sees that, at least in this calculation, the addition 
of processes which mix hard and soft partons do not destroy the agreement for the 
model with V2/n which is presumably a soft process. 

7.5 Hadron Formation Time 

In the discussion of the suppression of pions for px > 2 GeV/c, we treat the 
pions as resulting from the fragmentation of hard-scattered quarks and gluons. 
The explanation of this suppression in terms of partonic energy loss assumes 
that the hadronic wave function only becomes coherent outside the medium. 
Protons have a different hadronic structure and larger mass, and so may have 
a different, shorter time scale for coherence. 

Following [220], we can estimate the formation time for the different mass 
hadrons at moderate pt in two different ways. According to the uncertainty 
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principle, the formation time in the rest frame of the hadron can be related 
to the hadron size, Rh- In the laboratory frame, the hadron formation time is 
then given by 



Tf W R h 



E h 



(33) 



where Rh is taken to be 0.5-1 fm. For a 10 GeV/c pion, this gives a formation 
time of 35-70 fm/c. For the pt =2.5 GeV/c pions considered in this section, 
the formation time is 9-18 fm/c, well outside the collision region. However for 
Pt =2.5 GeV/c protons, the corresponding formation time is only 2.7 fm/c 
in the vacuum, suggesting the possibility that the hadronization process may 
begin inside the medium. However the formation of such heavy particles would 
presumably be delayed in a deconfined medium until the entire system began 
to hadronize. 

If quarks and antiquarks from gluon splitting are assumed to combine into 
dipole color singlets leading to the final hadrons, the formation time may be 
estimated from the gluon emission time. Then the formation time for a hadron 
carrying a fraction z of the parton energy is given by 



If z is 0.6-0.8 and kx ~ Aqcd, proton formation times in the range of 1-4 
fm result [220]. Such values again imply formation of the proton within the 
medium. Thus, it is possible that differing (and perhaps complicated) inter- 
actions with the medium may produce different scalings of proton and pion 
production and result in modified fragmentation functions in Au+Au colli- 
sions. However, most expectations are that this should lead to greater sup- 
pression rather than less. In fact, modified fragmentation functions measured 
in electron deep-inelastic scattering on nuclei by the HERMES experiment are 
often interpreted in terms of additional suppression for hadrons forming in the 
nuclear material. 

7.6 Hard-Scattering Physics 

If the dominant source of (anti)protons at intermediate pt is not soft physics, is 
the explanation a medium-modified hard-process source? The near-side part- 
ner yields indicate that a significant fraction of the baryons have jet-like part- 
ners. However, in the parton energy loss scenario as described in Section 6, 
hard-scattered partons lose energy in medium prior to hadronization. Thus 
one would expect the same suppression for baryons and mesons. Furthermore, 




(34) 
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we know that the (anti)protons £1X6 clS suppressed above pt = 5 GeV/c in a 
manner similar to pions . Hence for this explanation to be correct, there must 
exist a mechanism by which only partons leading to baryons between 2 and 5 
GeV/c in p T escape suppression. 

Another key piece of information is that the elliptic flow v 2 for protons is large 
for pt in the range 2-4 GeV/c. At low pt this collective motion is attributed to 
different pressure gradients along and perpendicular to the impact parameter 
direction in semi-central collisions. At higher px it has been hypothesized that 
one could observe a v 2 due to smaller partonic energy loss for partons traveling 
along the impact parameter direction (shorter path in the medium) as opposed 
to larger partonic energy loss in the perpendicular direction (larger path in the 
medium) . However, the data suggest that the pions have a large energy loss (a 
factor of five suppression in central Au+Au reactions) , while the protons do 
not. In this case one might expect that if the source of proton v 2 were energy 
loss, then proton v 2 would be significantly less than the v 2 for the pions. In 
fact, the opposite is experimentally observed: for px > 2 GeV/c, the proton 
v 2 is always larger than the pion v 2 . 

The contradictions the data create for both the "soft"- and "hard" -physics ex- 
planations may indicate that the correct physics involves an interplay between 
the two. 



7. 7 Conclusions 

The anomalous enhancement of (anti)protons relative to pions at intermedi- 
ate pt = 2-5 GeV remains a puzzle. At lower transverse momentum particle 
production is a long-wavelength "soft" process and the transport of these 
hadrons and their precursor partons is reasonably described by hydrodynam- 
ics. As observed at lower energies, soft particles emitted from an expanding 
system receive a collective velocity boost to higher pt resulting in an enhanced 
p/n and pjn ratio relative to proton-proton reactions at the same energy. We 
observe a similar phenomena at RHIC, for which the (anti) proton spectra and 
v 2 are roughly described in some hydrodynamic models up to approximately 
2 GeV/c. Another class of calculations, referred to as recombination models, 
also boosts soft physics to higher pt by coalescence of "dressed" partons. In 
the hydrodynamic models the quanta which are flowing are initially partons 
and then hadrons. The recombination models describe comoving valence par- 
tons which coalesce into hadrons, and do not reinteract. These two points of 
view may not be entirely contradictory, since both include a flowing partonic 
phase. If fact, it may be that the recombination models provide a mechanism 
by which hydrodynamics works to a much higher pt than one might expect. 
The simplifying assumption of hadrons which do not interact is most probably 
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an oversimplification and further refinement of the models will include this, 
though it may be that the hadronic phase will not modify the spectra as much 
as the hydrodynamic models might predict. 

In both models, the (anti)proton enhancement as a function of centrality can 
be tuned to reproduce the apparent binary collision scaling observed in the 
data. An important distinction between the two is that in one case this en- 
hancement is mass dependent and in the other it comes from the combination 
of quark momenta and thus distinguishes between baryons and mesons 10 . 
Rep for the is similar to other mesons despite the fact that they are more 
massive than protons. This scaling with quark content, as opposed to mass, 
favors recombination models. 

Further investigations into these intermediate px baryons reveals a near-angle 
correlation between particles, in a fashion characteristic of jet fragmentation. 
The near-angle associated particle yield increases by almost a factor of two in 
going from proton-proton and deuteron-gold reactions to gold-gold peripheral 
collisions. In addition, the partner yield is similar for trigger pions and protons, 
except in the most central gold-gold reactions. This appears to indicate a hard 
process source for a significant fraction of these baryons in contrast to the 
previous mentioned physics scenarios. Quantifying the precise contribution is 
an important goal for future measurements. 

The large (anti) baryon to pion excess relative to expectations from parton 
fragmentation functions at intermediate px — 2 — 5 GeV/c remains one of 
the most striking unpredicted experimental observations at RHIC. The data 
clearly indicate a new mechanism other than universal parton fragmentation as 
the dominant source of baryons and anti-baryons at intermediate pt in heavy 
ion collisions. The boosting of soft physics, that dominates hadron production 
at low px, to higher transverse momentum has been explored with the con- 
text of hydrodynamic and recombination models. However, investigations into 
these intermediate pt baryons reveals a near-angle correlation between parti- 
cles, in a fashion characteristic of jet fragmentation. If instead these baryons 
have a partonic hard scattering followed by fragmentation source, this frag- 
mentation process must be significantly modified. It is truly remarkable that 
these baryons have a large v 2 (typically 20%) indicative of strong collective 
motion and also a large "jet-like" near-side partner yield. At present, no the- 
oretical framework provides a complete understanding of hadron formation in 
the intermediate p T region. 



A caveat to this fact is that in the recombination models, it is the constituent- 
quark mass that is important, thereby giving a slightly larger mass to the strange 
quark. 
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8 FUTURE MEASUREMENTS 



The previous sections have documented the breadth and depth of the PHENIX 
data from the first three years of RHIC operations, along with the physics 
implications of those results. Here we describe those measurements required 
to further define and characterize the state of matter formed at RHIC. In 
particular, we note that the study of penetrating probes, which are the most 
sensitive tools in this endeavor, is just beginning. The PHENIX experiment 
was specifically designed to address these probes with capabilities that are 
unique within the RHIC program and unprecedented in the field of relativistic 
heavy ion physics. 

One can distinguish two broad classes of penetrating probes: 

(1) Hard probes created at the very early stage of the collision which propa- 
gate through, and could be modified by, the medium. These are the QCD 
hard-scattering probes and the main observables are high-p T particles 
coming from the fragmentation of jets, hidden charm (J/ip production), 
open charm and eventually also bottom quark and T production. 

(2) Electromagnetic probes (either real or virtual photons) which are created 
by the medium. Due to their large mean free path these probes can leave 
the medium without final-state interaction thus carrying direct informa- 
tion about the medium's conditions and properties. The main observables 
here are low-mass e + e~ pairs and the thermal radiation of the medium. 

By their very nature, penetrating probes are also rare probes and consequently 
depend on the development of large values of the integrated luminosity. In 
the present data set the reach for high-p^ particles in PHENIX extends to 
roughly 10 GeV/c, and lower-cross-section measurements such as charmonium 
are severely limited. The dramatic improvement of the machine performance 
in the year 2004 run provides confidence that both this data set and those 
from future RHIC runs will dramatically extend our reach in the rare probes 
sector. 

As part of a decadal planning of the RHIC operation, PHENIX has prepared 
a comprehensive document that outlines in great detail its scientific goals and 
priorities for the next 10 years together with the associated detector upgrade 
program needed to achieve them. The decadal plan [243] is centered around 
the systematic study of the penetrating probes listed above. The program 
is broad and can accommodate additions or modifications provided that a 
compelling physics case can be made. Measurements are mainly planned in 
Au+Au collisions at the full RHIC energy but they will be supplemented by 
other measurements varying the energy and/or the species and by the neces- 
sary reference measurements of p + p and p + A collisions. A short summary 
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is given below. 



8. 1 High-pT Suppression and Jet Physics 



The most exciting results to date at RHIC are the discovery of high-p^ suppres- 
sion of mesons, interpreted in terms of energy loss of quarks in a high-density 
medium, and the nonsuppression of baryons or equivalently, the anomalously 
high p/n ratio which still awaits a clear explanation. These two topics were 
extensively discussed in Sections 6 and 7, respectively. 

The data collected so far are superb. However, they suffer from limited reach in 
transverse momentum, limited particle identification capabilities and limited 
statistics in particular for detailed studies of jet correlations. PHENIX has 
a program for further studies of the high-p^-suppression phenomena and jet 
physics which aims at overcoming these limitations. 

It will be necessary to trace the suppression pattern to much higher p T to 
determine whether (and if so, when) the suppression disappears and normal 
perturbative QCD behavior sets in. High-luminosity runs will be needed, with 
at least a factor of 50 more statistics. PHENIX is particularly able to per- 
form these measurements with its excellent capability of triggering on high- 
momentum 7r°'s. 

PHENIX has performed several particle correlation analyses and has demon- 
strated that the experiment's aperture at mid-rapidity is sufficient to conduct 
these studies. Currently, these analyses are limited by the available statistics. 
Again, increasing the data sample by a factor of 50-100 will allow a variety of 
correlation studies using trigger particles with much-higher-momentum than 
studied to date. A particularly interesting case is the study of high-momentum 
7-jet correlations, which have vastly reduced trigger bias, since the trigger pho- 
tons propagate through the medium with a very long mean free path. 

To further elucidate the baryon puzzle, additional data is required with better 
separation between baryons and mesons. An upgrade consisting of an aerogel 
Cerenkov counter and a high-resolution TOF detector is expected to be com- 
pleted in time for the year 2006. A portion of this aerogel counter was already 
installed prior of the year 2004 run and performed according to expectations. 
Once completed, this high-p^ detector will allow identification of n,K/p to 
beyond 8 GeV/c in p T . 
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8.2 J/ip Production 



Suppression of heavy quarkonia is one of the earliest and most striking pro- 
posed signatures of deconfmement. The suppression mechanism follows di- 
rectly from the Debye screening expected in the medium, which reduces the 
range of the potential between charm quark and anti-quark pairs [244]. The 
NA38 and NA50 experiments have carried out a systematic study of J / ip and 
ip' at the CERN-SPS in p + p, p + A, light ion, and Pb+Pb collisions provid- 
ing some of the most intriguing results of the relativistic heavy ion program 
for more than ten years. The NA50 experiment observed an anomalous sup- 
pression of J/ip in central Pb+Pb collisions at y/s NN = 17.2 GeV [245]. The 
suppression, which is of the order of 25% with respect to the normal suppres- 
sion in nuclear matter, has been interpreted by the NA50 authors as evidence 
for deconfmement of quarks and gluons. Although this interpretation is not 
universally shared [246,247], the results of NA38 and NA50 demonstrate the 
utility and great interest in understanding the fate of charmonium in dense 
nuclear matter. 



The theoretical expectations at RHIC energies are not at all clear. They range 
from total suppression in the traditional Debye screening scenario to enhance- 
ment in coalescence models [248,249,250] and in statistical hadronization mod- 
els [251,252], of c and c quarks. Although some versions of the coalescence 
model seem disfavored from our very limited data set [51], a more conclusive 
statement on these models has to await the much larger data set of the year 
2004 run. 



PHENIX has unprecedented capabilities for the study of the J/ip in Au+Au 
collisions. The J/ip can be measured via its decay channel at forward 

and backward rapidities in the muon spectrometers and via its e + e~ decay 
channel at mid-rapidity in the central arm spectrometers. From the recorded 
luminosity of the year 2004 run, we expect several thousand and ~500 J / in 
the muon and central arms, respectively. This data set will allow us a first look 
at the J/ip production pattern at RHIC. However, it could well be marginal 
for a complete characterization as a function of centrality and px, so that it is 
likely that further higher-luminosity runs will be required. Also the p + p and 
d+Au baseline measurements performed in the year 2001 - 2003 runs have 
large statistical uncertainties, and higher-statistics versions for these colliding 
species will be needed. A high-luminosity p+p run is planned in the year 2005 
and high-luminosity d + A or p + A are still to be scheduled in the next years. 
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8.3 Charm Production 



Charm quarks are expected to be produced in the initial hard collisions be- 
tween the incoming partons. The dominant mechanism is gluon fusion and 
thus the production cross section is sensitive to the gluon density in the initial 
state. The cc production cross section is sizable at RHIC energies with a few cc 
pairs and therefore several open charm mesons per unit of rapidity in central 
Au+Au collisions. As a result, charm observables become readily accessible at 
RHIC and offer additional and extraordinarily valuable diagnostic tools. For 
example, it is vitally important to perform measurements of charm flow and 
to determine the energy loss of charm quarks in the medium. Such measure- 
ments will determine if the bulk dynamics observed for light quarks extend to 
charm quarks, which could in fact have very different behavior due to their 
much larger mass. Again the potential of PHENIX in this domain is unique 
with its capability of measuring open charm in a broad rapidity range, in the 
central and muon arms, via both the electron and muon decay channels. An 
additional unique feature is the possibility to measure correlated semileptonic 
charm decays by detecting e — \x coincidences from correlated DD decays. 
Such a measurement is particularly interesting for the study of charm-quark 
energy loss which may differ significantly from that observed for lighter quarks 
[253,254,255]. A first study of e — /j, coincidences should be feasible with the 
year 2004 data. 

To date PHENIX has measured charm production cross section in an indirect 
way through high-p T single electrons [42,57] assuming that all electrons (after 
measuring and subtracting the contributions from light hadrons and photon 
conversions) originate from the semileptonic decays of charm quarks. Although 
the charm cross section has large uncertainties, the centrality dependence of 
the charm rapidity density demonstrates that charm production follows binary 
scaling as shown in Fig. 34. Improvements and additional information are 
expected from the much higher statistic of the year 2004 data. 

A qualitatively new advance for PHENIX in the charm and also the beauty 
sector will be provided by the implementation of the silicon vertex detector. 
An upgrade project is underway to install in the next five years a silicon 
vertex tracker, including a central arm barrel and two end caps in front of the 
two muon spectrometers. The vertex tracker will allow us to resolve displaced 
vertices and therefore to directly identify open charm mesons via hadronic, 
e.g. D — > Kit, as well as semi-leptonic decays. The heavy-quark physics topics 
accessible with the vertex tracker include production cross section and energy 
loss of open charm and open beauty, and spectroscopy of charmonium and 
bottomonium states, each of which should provide incisive new details on the 
properties of the created medium. 
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8.4 Low-Mass Dileptons 



Low-mass dileptons are considered the most sensitive probe of chiral symmetry- 
restoration primarily through p meson decays. Due to its very short lifetime 
(r = 1.3 fm/c) compared to that of the typical fireball of ~ 10 fm/c, most 
of the p mesons decay inside the medium providing an unique tool to observe 
in-medium modifications of its properties (mass and/or width) which could be 
linked to chiral symmetry restoration. The situation is somewhat different but 
still interesting for the uj and mesons. Because of their much longer lifetimes 
(r = 23 fm/c and 46 fm/c for the uo and 0, respectively ) they predominantly 
decay outside the medium, after regaining their vacuum properties, with only 
a small fraction decaying inside the medium. Since the measurement integrates 
over the history of the collision, this may result in a small modification of the 
line shape of these two mesons which PHENIX might be able to observe with 
its excellent mass resolution. PHENIX also has the unprecedented capability 
of simultaneously measuring within the same apparatus the meson decay 
through e + e~ and K + K~ channels. The comparison of the branching ratios to 
these two channels provides a very sensitive tool for in-medium modifications 
of the and K mesons. 

The CERES experiment at CERN has confirmed the unique physics potential 
of low-mass dileptons [256,257,258]. An enhancement of electron pairs was ob- 
served in the mass region m = 0.2-0.6 GeV/c 2 in Pb+Au collisions at ^snn = 
17.2 GeV with respect to p+p collisions. The results have triggered a wealth of 
theoretical activity and can be explained by models which invoke in-medium 
modification of the p meson (dropping of its mass and/or broadening of its 
width) [259]. The precision of the CERES data has been so far insufficient to 
distinguish between the different models. Results with higher statistics and 
better mass resolution are expected from the NA60 experiment that is study- 
ing the production of low- mass dimuons in In+In collisions [260]. Theoretical 
calculations [261] show that the enhancement should persist at RHIC energies 
and that PHENIX with its excellent mass resolution has an unique opportu- 
nity to do precise spectroscopy of the light vector mesons and to shed more 
light on the origin of the enhancement of the low-mass-pair continuum. 

The measurement of low-mass electron pairs is however a very challenging 
one. The main difficulty stems from the huge combinatorial background cre- 
ated by the pairing of e + and e~ tracks from unrecognized 7r° Dalitz decays 
and 7 conversions. PHENIX is developing a novel Cerenkov detector that, in 
combination with the recently installed coil which makes the magnetic field 
zero close to the beam axis, will effectively reduce this combinatorial back- 
ground by almost two orders of magnitude [262]. The detector, operated in 
pure CF 4 , consists of a 50-cm-long radiator directly coupled, in a windowless 
configuration, to a triple GEM detector which has a Csl photocathode evap- 
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orated on the top face of the first GEM foil and pad read out at the bottom 
of the GEM stack [263]. The R&D phase to demonstrate the validity of the 
concept is nearing completion. The detector construction phase is starting 
now with installation foreseen in time for the year 2006 — 2007. With this 
detector PHENIX will have the unprecedented ability to perform high-quality 
measurements over the whole dilepton mass range from the n° Dalitz decay 
up to the charmonium states. 



8.5 Thermal Radiation 

A prominent topic of interest in the field of relativistic heavy-ion collisions 
is the identification of the thermal radiation emitted by the system and in 
particular the thermal radiation emitted by the quark-gluon plasma via qq 
annihilation. Such radiation is a direct fingerprint of the matter formed and 
is regarded as a very strong signal of deconfinement. Its spectral shape should 
provide a direct measurement of the plasma temperature. 

In principle the thermal radiation can be studied through real photons or 
dileptons, since real and virtual photons carry basically the same physics mes- 
sage. In practice the measurements are extremely challenging. The thermal 
radiation is expected to be a small signal compared to the large background 
from competing processes, hadron decays for real photons and Dalitz decays 
and 7 conversions for dileptons, the former being larger by orders of mag- 
nitude compared to the latter. But in both cases, a very precise knowledge 
of all these sources is an absolutely necessary prerequisite. After subtracting 
these sources, one still needs to disentangle other contributions which might be 
comparable or even stronger, mainly the contributions of initial hard-parton 
scattering to direct photons and of semileptonic decays of charm mesons to 
dileptons. 

Theoretical calculations have singled out the dilepton mass range m — 1- 
3 GcV/c 2 as the most appropriate window where the QGP radiation could 
dominate over other contributions [264,265]. Measurements in this interme- 
diate mass range carried out at the CERN SPS by HELIOS and NA50 have 
revealed an excess of dileptons, but this excess could be explained by hadronic 
contributions [266]. 

There is no conclusive evidence for QGP thermal photons from the CERN 
experiments (for a recent review see [267]). From the theoretical point of view 
it is clear that in the low-p-r region (p? < 2 GeV/c) the real photon spectrum 
is dominated by hadronic sources and the thermal radiation from the hadron 
gas. It is only in the high-p-^ region where one might have a chance to observe 
the thermal radiation from the QGP. 
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PHENIX has measured direct real photons at pt > 4 GeV/c from the initial 
hard scatterings [73]. The errors are relatively large leaving room for a compa- 
rable contribution of thermal photons. The high statistics of the year 2004 run 
will provide the first real opportunity to search for the QGP thermal radia- 
tion in PHENIX both in the dilepton and real photon channels. However, the 
search for this elusive signal might take some time as it will probably require 
equally-high-statistics runs of reference data in p + p and p + A collisions for 
a precise mapping of all the other contributions (hadronic + pQCD for real 
photons and hadronic + charm for dileptons). 



9 SUMMARY AND CONCLUSIONS 



The PHENIX data set from the first three years of RHIC operation provides an 
extensive set of measurements, from global variables to hadron spectra to high- 
er physics to heavy-flavor production. From this rich menu we have reviewed 
those aspects of the present data that address the broad features of the matter 
created in Au+Au collisions at RHIC, namely, energy and number density, 
thermalization, critical behavior, hadronization, and possible deconfinement. 

We first investigated whether the transverse energy and multiplicity measure- 
ments of PHENIX demonstrate that a state of high-energy-density matter is 
formed in Au+Au collision at RHIC. We estimated from our dE T /drj measure- 
ment that the peak energy density in the form of created secondary particles 
is at least 15 GeV/fm 3 . If we use a thermalization time of 1 fm/c provided by 
the hydrodynamic models from the elliptic flow, then the value of the energy 
density of the first thermalized state would be in excess of 5 GeV/fm 3 . These 
values are well in excess of the ~1 GeV/fm 3 obtained in lattice QCD as the 
energy density needed to form a deconfined phase. Naive expectations prior 
to RHIC turn-on that dE T /drj and dN ch /dr] could be factorized into a "soft" 
and a pQCD jet component are not supported by the data. Results from a 
new class of models featuring initial-state gluon saturation compare well with 
RHIC multiplicity and E T data. 

We then examined our data and various theoretical models to investigate the 
degree to which the matter formed at RHIC appears to be thermalized. The 
measured yields and spectra of hadrons are consistent with thermal emission 
from a strongly expanding source, and the observed strangeness production 
is consistent with predictions based on complete chemical equilibrium. The 
scaling of the strength of the elliptic flow v 2 with eccentricity shows that a 
high degree of collectivity is built up at a very early stage of the collision. 
The hydro models which include both hadronic and QGP phases reproduce 
the qualitative features of the measured ^(pr) of pions, kaons, and protons. 
These hydro models require early thermalization (r therm < 1 fm/c) and high 
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initial energy density e > 10 GeV/fm 3 . These points of agreement between the 
data and the hydrodynamic and thermal models can be interpreted as strong 
evidence for formation of high-density matter that thermalizes very rapidly. 

However several of the hydro models fail to reproduce the ^(pt) of pions, pro- 
tons, and spectra of pions and protons simultaneously. Given this disagreement 
it is not yet possible to make an unequivocal statement regarding the pres- 
ence of a QGP phase based on comparisons to hydrodynamic calculations. 
The experimentally measured HBT source parameters, especially the small 
value of -Riong and the ratio -R ou t/-Rsidc ~ 1, are not reproduced by the hy- 
drodynamic calculations. Hence we currently do not have a consistent picture 
of the space-time dynamics of reactions at RHIC as revealed by spectra, t>2 
and HBT. These inconsistencies prevent us from drawing firm conclusions on 
properties of the matter such as the equation of state and the presence of a 
mixed phase. 

Critical behavior near the phase boundary can produce nonstatistical fluc- 
tuations in observables such as the net-charge distribution and the average 
transverse momentum. Our search for charge fluctuations has ruled out the 
most naive model of charge fluctuations in a QGP, but it is unclear if the 
charge fluctuation signature can survive hadronization. Our measurement of 
(Pt) fluctuations is consistent with the effect expected of high-pr jets, and it 
gives a severe constraint on the fluctuations that were expected for a sharp 
phase transition. 

Many of these observables — for instance, large dE/drj and dN c h/dr], strangeness 
enhancement, strong radial flow, and elliptic flow — have been observed in 
heavy ion collisions at lower energies. We have found smooth changes in these 
observables as a function of ^snn from AGS energies to SPS energies to RHIC 
energy. The dE T /drj increases by about 100% and the strength of the elliptic 
flow increase by about 50% from SPS to RHIC. The strangeness suppression 
factor 7 S and the radial expansion velocity {(3t) vary smoothly from AGS to 
RHIC energies. No sudden change with collision energy has been observed. 

The strong suppression of high-p T particle production at RHIC is a unique phe- 
nomenon that has not been previously observed. Measurements of two-hadron 
azimuthal-angle correlations at high px and the xt scaling in Au+Au collisions 
confirm the dominant role of hard scattering and subsequent jet fragmenta- 
tion in the production of high-p^ hadrons. Measurements in deuteron-gold 
collisions demonstrate that any initial-state modification of nuclear parton 
distributions causes little or no suppression of hadron production for p T > 2 
GeV/c at mid-rapidity. This conclusion is further strengthened by the observed 
binary scaling of direct photon and open charm yields in Au+Au. Combined 
together, these observations provide direct evidence that Au+Au collisions at 
RHIC have produced matter at extreme densities. 
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Medium- induced energy loss, predominantly via gluon bremsstrahlung emis- 
sion, is the only currently known physical mechanism that can fully explain 
the magnitude of the observed high-p^ suppression. The approximately flat 
suppression factor Raa{pt) observed in the data, which was predicted by the 
GLV energy loss model, rules out the simplest energy loss models which pre- 
dicted a constant energy loss per unit length. However, the model by Wang 
et al. obtains the same flat Raa{pt) from apparently different physics. From 
the GLV model, the initial gluon number density, dn g /dy ~ 1000 and initial 
energy density, e 15GeV/fm 3 , have been obtained. These values are consis- 
tent with the energy density obtained from our dE T /drj measurement as well 
as ones from the hydro models. 

The large (anti)baryon to pion excess relative to expectations from parton 
fragmentation functions at intermediate pt (2 — 5 GeV/c) is both an un- 
predicted and one of the most striking experimental observation at RHIC. 
The data clearly indicates that a mechanism other than universal parton frag- 
mentation is the dominant source of (anti-)baryons in the intermediate pt 
range in heavy ion collisions. The boosting of soft physics to higher trans- 
verse momentum has been explored within the context of hydrodynamics and 
recombination models. Hydrodynamic models can readily explain the baryon 
to meson ratio as a consequence of strong radial flow, but these models have 
difficulties reproducing the difference in v<i between protons and mesons above 
2 GeV/c. Recombination models provide a natural explanation for the large 
baryon to meson ratio as well as the apparent quark-number scaling of the 
elliptic flow. However, investigations into these intermediate pt baryons re- 
veal a near-angle correlation between particles, in a fashion characteristic of 
jet fragmentation. If instead these baryons have a partonic hard scattering 
followed by fragmentation, this fragmentation process must be significantly 
modified. It is truly remarkable that these baryons have a large i>2 of ~ 20 % 
typically indicative of strong collective motion and also a large jet-like near- 
side partner yield. At present, no model provides a complete understanding of 
hadron formation in the intermediate px regime. 

The initial operation of RHIC has produced the impressive quantity of signifi- 
cant results described above. These striking findings call for additional efforts 
to define, clarify and characterize the state of matter formed at RHIC. Further 
study of the collisions using hard probes such as high-p^ particles, open charm, 
and J/ip, and electromagnetic probes such as direct photons, thermal photons, 
thermal dileptons, and low-mass lepton pairs are particularly important. The 
utilization of these penetrating probes is just beginning, and we expect these 
crucial measurements based on the very-high-statistics data of the year 2004 
run will provide essential results towards understanding of the dense matter 
created at RHIC. 

Advances in the theoretical understanding of relativistic heavy ion collisions 
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is vital for the quantitative study of the dense matter formed at RHIC. While 
there is rapid and significant progress in this coherent and consistent 

picture of heavy ion collisions at RHIC, from the initial formation of the 
dense matter to the thermalization of the system to the hadronization to 
the freezeout, remains elusive. With such a consistent model, it will become 
possible to draw definitive conclusions on the nature of the matter and to 
quantitatively determine its properties. The comprehensive data sets from 
global variables to penetrating probes provided by PHENIX at present and 
in the future will prove essential in constructing and constraining a consistent 
model of heavy ion collisions to determine the precise nature of the matter 
created at RHIC. 

In conclusion, there is compelling experimental evidence that heavy-ion col- 
lisions at RHIC produce a state of matter characterized by very high energy 
densities, density of unscreened color charges ten times that of a nucleon, 
large cross sections for the interaction between strongly interacting particles, 
strong collective flow, and early thermalization. Measurements indicate that 
this matter modifies jet fragmentation and has opacity that is too large to 
be explained by any known hadronic processes. This state of matter is not 
describable in terms of ordinary color-neutral hadrons, because there is no 
known self-consistent theory of matter composed of ordinary hadrons at the 
measured densities. The most economical description is in terms of the under- 
lying quark and gluon degrees of freedom. Models taking this approach have 
scored impressive successes in explaining many, but not all, of the striking fea- 
tures measured to date. There is not yet irrefutable evidence that this state of 
matter is characterized by quark deconfinement or chiral symmetry restora- 
tion, which would be a direct indication of quark-gluon plasma formation. The 
anticipated program of additional incisive experimental measurements com- 
bined with continued refinement of the theoretical description is needed to 
achieve a complete understanding of the state of matter created at RHIC. 
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